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VTOL CONTROLS FOR SHIPBOARD LANDING 
by 

Christopher Graham McMuldroch 


Submitted to the Department of Aeronautics and Astronautics on May 11, 
1979, in partial fulfillment of the requirements for the degree of Master 
of Science. 


ABSTRACT 

The problem of landing a VTOL aircraft on a small ship in rough 
seas using an automatic controller is examined. The controller design 
uses the linear quadratic gaussian results of modern control theory. 
Performance bounds and control requirements are obtained for two control 
concepts . 

V 

Linear time invariant dynamic models are developed for the aircraft, 
ship and wave motions. The VTOL is modelled at stationary hover, sup- 
ported by a lift fan in the nose and a pair of lift/cruise fans at the 
wing roots. Thrust control dynamics are modelled as first order, while 
aerodynamic effects are omitted. The ship is modelled at cruising speed 
in waves typical of sea state five conditions. The wave model, which 
drives ship motion, is shaped white noise with a theoretical Neumann Spec- 
trum. All states are assumed measurable without reconstruction or noise. 

The first control concept is a "hover controller", which commands 
the aircraft to track position and orientation of the ship deck. An 
important goal is to use only low levels of control power. Commands for 
this task are generated by the solution of the steady state linear quad- 
ratic gaussian regulator problem. Analytical performance and control 
requirement tradeoffs are obtained. 

The second control concept is a "landing controller", which is 
designed to command the aircraft from stationary hover along a smooth, 
low control effort trajectory, to a touchdown on a predicted crest of 
ship motion. The design problem in this case is formulated and solved 
as an approximate finite-time linear quadratic stochastic regulator. 
Performance and control results are found by Monte Carlo simulations. 


Thesis Supervisor: Gunter Stein, Ph.D. 

Title: Adjunct Professor of Electrical 

Engineering 
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- ocean wave model driving noise vector 

- low order ship model input for pitch mode 

- low order ship model input for roll mode 

- noise input vector 

- input noise vector for low order ship model 

- standard deviation or 

- width constant in 

- standard deviation of v 

- rms value of z. _ 
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- time variable 

- aircraft roll angle 

- ship roll angle 

- aircraft yaw angle 

- ship yaw angle 

- angular frequency 

- frequency of encounter of ocean waves by moving ship 

- ocean wave model ith natural frequency 

- natural heave frequency of ship 
-ith natural frequency in wave model or 

- ith natural frequency in low order ship model transfer 
function for mode i 

- natural pitch frequency of ship 

- natural roll frequency of ship 

- frequency of ocean waves seen by stationary observer 

- aircraft angular velocity in the aircraft body frame 

- derivative with respect to time of (•) 

] partial derivative of (•) with respect to [•] 

- transpose of («) 

- time derivative of (•) 

- absolute value of (*) 


- percent 
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CHAPTER 1 
INTRODUCTION 

In recent years, interest in vertical take-off and landing (VTOL) 
aircraft has been demonstrated by a number of studies on potential uses, 
designs, and technologies required for such vehicles. 

Among their many uses, one of the most interesting and significant is 
the operation of VTOL vehicles from ships at sea. These aircraft could be 
carried on smaller ships possessing far fewer resources than conventional 
aircraft carriers. By virtue of combining the hover capability of rotary 
wing aircraft and the cruise efficiency of fixed wing vehicles, VTOL air- 
craft would bring the benefits of long range air support to ships previously 
only capable of carrying helicopters. 

Several VTOL aircraft design concepts have been proposed (1,2, 3, 4, 5). 
One of these is a small lift/cruise fan vehicle supported in hover by 
three thrust fans: one in the aircraft nose, and the other two at the wing 
roots. These fans are driven by three gas turbine engines interconnected 
by shafts or gas ducts in a manner permitting power exchange among the fans 
without increasing engine power. This provides engine-out safety and rapid 
vehicle attitude control. 

One of the critical problems of ship-borne VTOL operation is the 
control of landing under severe wave and wind conditions. The vehicles 
must be capable of achieving soft touchdowns with very l imi ted control 
authority (thrust to weight capability) under motion conditions as 
severe as "sea state five" (.12 ). Roughly speaking, this corresponds 
to wave heights of six feet or more and wind speeds as high as twenty 
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knots. This thesis investigates two concepts for automatic control of 
landings under these conditions. 

1.1 Background References and Previous Work 

A qualitative introduction to recent VTOL ideas is found in refer- 
ences (1, 2, 3, 4, 5). The first four of these describe various V/STOL 
vehicles and concepts. The last also describes V/STOL requirements 
including shipboard landings. More detailed papers on a wide range of 
V/STOL topics are located in reference ( 6 ) . 

The results of a landing controller design and evaluation, which 
include a manually piloted shipboard landing task, are reported in 
reference ( 7 ), "Design and Evaluation of an Integrated Flight Control 
System Concept for Manual IFR VTOL Operations", by Vernon K. Merrick 
and Ronald M. Gerdes. The basic shipboard landing problem in this 
reference is similar to the one investigated here. In both cases, the 
goal is to land VTOL aircraft accurately and gently aboard a ship at 
sea. However, there are some notable differences-. First, the reference 
addresses a manually piloted landing problem, rather than an automatic 
one as done here. Secondly, the aircraft modelled was a transport 
aircraft with six thrust fans, not a small vehicle with only three fans 
as is used in this report. Next, the ship motions modelled in the 
reference are generated by a combination of sinusoidal components; 
whereas the ones modelled here are stochastic. Finally, the control 
design techniques in the reference were classical in contrast to the 
modern approach followed in this work. 
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1.2 Objective, Approach and Contributions 

\ 

The objective of this work is to design a landing controller to 
guide a VTOL aircraft from stationary hover along a smooth trajectory 
to a gentle landing on the ship deck. At the same time, the control 
authority used to follow the trajectory is to be kept small. The in- 
tent is to achieve these objectives with design techniques based on 
modern linear quadratic control. Besides outlining design techniques, 
the results are also intended to provide design guidelines and per- 
formance bounds for eventual practical landing control implementations. 

It is not a goal of this research to form a detailed engineering design. 

Three steps were used to reach these objectives. The first step 
was to construct mathematical models for each of the system's components, 
including VTOL aircraft, ship, and ocean waves. The second step was the 
design of a ship deck chasing "hover controller". The performance of 
this design motivated the third step, which was the design of a "landing 
controller". This second controller guides the aircraft to a predicted 
peak of ship deck motion. 

The aircraft modelled is a small VTOL with three lift fans whose 
thrust magnitude and direction can be varied to control translational 
and rotational motion. Ship motion is represented by three parallel 
second order systems for pitch, heave and roll, respectively. These 
three systems are driven by a pair of sixth order wave models, which 
in turn are driven by independent white noises . The wave models approxi- 
mate the theoretical Newmann wave spectrum. Ship model parameters are 
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chosen to represent a small ship of about 400 foot length in sea state 
five conditions. To date, pitch-heave dynamics only have been used in 
the control design studies. 

The hover control goal is to make the aircraft track the ship 
landing deck motion at constant vertical separation. Implemented for 
automatic landings, the vertical separation command would be replaced 
by a smoothly decreasing function. The controller design is a steady 
state linear quadratic regulator chosen to minimize weighted separation 
error and vertical thrust control effort. Curves of separation error 
versus vertical thrust control are generated. In addition, a curve for 
the required vertical sensor and control bandwidths versus control 
effort is determined. 

The thrust to weight requirements of the above controller indicate 
that an alternative approach is desirable for low control power landings. 
The approach taken is to predict ship deck motion forward to a peak. 

Then a smooth nominal trajectory is computed and followed to the pre- 
dicted peak point and time. The nominal trajectories are generated by 
a finite end-time linear quadratic regulator using a low order aircraft 
model. Monte Carlo simulations of this nonlinear landing algorithm 
show that the required control authority is reduced by a factor of 
almost two over that required by the hover controller. This demonstrates 
the potential utility of ship motion prediction. 

This study demonstrates how modern control theory can be applied 
to the VTOL shipboard landing problem. "First iteration" control designs 
have been generated which provide basic performance trade-off information 
for shipboard VTOL operation. Low touchdown velocities can be achieved 
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provided that the control bandwidth exceeds the ship motion spectrum 
peak at approximately one radian/ s econd . Moreover, ship motion pre- 
diction appears to offer a promising way to reduce the control require- 
ments. These findings are tentative, of course, because idealistic 
modelling and sensing assumptions were made in the design. However, 
the work has been structured so that more detailed investigations can 
be made of performance variations with aerodynamic effects, ground 
effects, sea state conditions, sensor constraints, human piloting 
considerations, and inclusion of full longitudinal and lateral dynamics. 

1.3 Overview of the Chapters 

The material in Chapters Two through Five corresponds directly 
with the synopsis presented above. Chapter Two presents the model. 
Linear, time- invariant aircraft, wave, and ship models are developed 
from the nonlinear equations and frequency domain information available 
in the literature. Chapter Three presents the hover controller design 
and evaluation. This controller is studied with the initial intention 
of generalizing it for landing control. The high control authority 
required there leads to the landing controller in Chapter Four. This 
more promising controller is developed and then evaluated by Monte 
Carlo simulations. Chapter Five presents general conclusions, and 
suggests avenues for further investigation. 
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CHAPTER 2 
MODELLING 

2.1 Introduction 

The physical system to be modelled has three subsystems: aircraft, 

ship, and waves. For each of these, a linearized model is derived based 
on a nominal condition and simplifying assumptions. This chapter dis- 
cusses the coordinate systems employed and the nominal configuration of 
aircraft, ship, and waves. Next, the representation of the overall model 
as the augmentation of its subsystems is presented. Following this 
general discussion, the model of each subsystem is derived. In Chapter 3, 
a low order ship motion model is used for comparison of results with the 
normal wave-ship model. This low order model is also derived in this 
chapter. 

2.2 Coordinate Systems 

Figures (2.2.1, 2, and 3) illustrate the coordinate systems used. 

The aircraft and ship each have two associated coordinate systems. Each 
of these are standard, orthogonal, right-handed reference frames. They 
are referred to as "body frames" fixed to the vehicle axes of symmetry 
and "earth frames" fixed to the earth vertical and horizontal directions . 
The x-axes all point nominally forward, the z-axes point nominally down- 
ward, and- the y-axes are oriented to the right completing the right- 
handed frames. Note the effect of the z-axes definition. Heights above 
the center of rotation have positive altitude and heave, but negative 
z-coordinate. For the aircraft, the body fixed coordinate system is 
centered at the center of mass. For the ship, the corresponding frame 


t 


34 



Side View of 

Figure 2. 2>1 
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Plan View of Aircraft-Ship Configuration 
Figure 2.2.2 
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frame is located at a point approximating the center df rotation. The 
earth reference frame for the ship moves with the ship center of rota- 
tion, but maintains a constant height and a steady orientation with the 
Xgg and y ES ~axes locally horizontal, z-axis vertical, and x-axis parallel 
to the ship velocity vector. The aircraft earth frame is fixed with 
respect to the ship-earth frame. It is centered at some constant desired 
height above the nominal position of the ship landing deck. 

The following is a summary of the axes labels and other geometry: 
Ship body axes *j S , / BS > * Bs i 

Ship earth axes x^, y ES , z ES ; 

Ship body heave axis 
Ship earth heave axis 

Ship landing point (SLP) is at ship coordinates 

<X SLP- W W “ here ^SLP = ° ! 

Ship pitch angle 6 ; 

J 

Ship roll angle <j> ; 

o 

Ship yaw angle $ = 0); 

Aircraft body area . y^, z^s 

Aircraft earth axes x^, y^, z^; 

Aircraft pitch angle 0 ; 

A 

Aircraft roll angle <p ; 

A 

Aircraft yaw angle \p ; 

A 

Aircraft position relative to the ship deck in earth 
coordinates z ; 

A* 3 

Aircraft pitch orientation relative to the ship 0 ; 

A- S 
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Heading angle of the wind over the ship deck 0 ; 

AW 

Heading angle of wind and waves over ocean 0^. 

The three drawings above and the notation just introduced are useful in 

describing the relative position of the vehicles. The ship is steaming 

into the wind and waves at an angle 8^, so that the wind over the deck 

crosses at an angle 0, r , (8.,, = 0 in the design work to date). This is 

AW Aw 

the heading angle maintained by the hovering aircraft in order to avoid 
excessive sideforces. The aircraft position, z^ g, and orientation, 0^_g, 
relative to the ship landing point (or pad) are quantities to be controlled. 

2.3 Overall Model 

2.3.1 Augmented Model Dynamical Equation 

The overall linearized model of the aircraft, ship and waves is 
written as an augmented matrix differential equation. This equation is 
given in Table (2. 3. 1.1) to indicate the augmented nature of the model. 

The notation used for this equation is explained below: 

x, x are the overall model state vector and its derivative; 

_u is the vector of aircraft controls; 

is the vector of noises driving the wave model; 

A is the linearized, open loop, augmented or overall system matrix; 

B is the overall system control input matrix; 

is the overall system plant noise input matrix; 
x^, x< 3 > -^OWPH’ -^OWR are t * ie subsystem state vectors for the air- 
craft, ship, ocean wave pitch-heave model, and ocean wave roll 



39 


THE OVERALL LINEARIZED OPEN LOOP AUGMENTED MODEL EQUATION 


• X '• = ■ 

Ax + 

Bu + 

B N £ OW 

A 1 
a a 

0 

0 

1 0 

o 1 

A 

A s 

A 

SWPH 

1- - A - - - 
a swr 

0 1 

0 

A ow 

1 0 


0 

0 

1 A ow 


b a 


0 

0 

0 


0 

0 

0 

b n 3 

1 1 

is 

j 

i i 

0 

0 


O 

B ow 


-A 



1 

1 

tol 

*1 

1 


P OWPH 

_ gpWPH 

^OW 

P OWR 

— OWR 




TABLE 2. 3. 1.1 
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model, respectively; 

p , p are the pitch-heave and roll ocean wave model noise 
UWK 

inputs ; 


A , A , A are the aircraft, ship and the ocean wave model system 

A. j UW 

matrices: note pitch-heave and roll waves have the same system 

matrix; 


B , B are the aircraft and ocean wave model input matrices. 

A Vj W 

Again B QW is shared by the two models; 


^SWPH’ ^SWR are cou P^ n § matrices driving ship motion from the 
pitch-heave and roll wave models respectively. 


The elements of the subsystem matrices are determined in the sections 
that follow. In the standard model with the standard wave-ship models, 
there are a total of 36 states. The equation for the overall model 
using the low order ship model is shown in Table (2. 3. 1.2). This model 
is driven by p^, which acts directly on the ship model. The low order 


model has 24 states. 



41 


THE OVERALL LINEARIZED OPEN LOOP MODEL USING THE LOW ORDER SHIP MODEL 
x. + Ax + Bu + 



TABLE 2. 3. 1.2 
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2.3.2 Overall Model Output Equation 


The overall model output is basically the augmentation of all the 
individual model outputs; however, in addition, other outputs are included. 
From the aircraft and ship outputs of position and orientation are formed 
the additional outputs of relative aircraft-ship vertical separation and 
pitch orientation, as well as their relative rates. This is done by 
differencing the appropriate outputs of the aircraft and ship. In Table 
(2. 3. 2.1), these additional outputs are denoted ant * are obtained by: 


*A-S 



( 2 . 3 . 2 . 1 ) 


In the table, the following notation is used in addition to that of the 
previous section: 

£ is the overall model output vector; 

C is the output matrix; 

X&’ .Zs* 2 owPH’ -^OWR are the ouCput vectors for the aircraft, 
ship, pitch-heave, and roll ocean wave models, respectively. 

C A’ C S’ C 0W are the out P ut matrices of the aircraft, ship, 
and two identical ocean wave models . 

This completes the discussion of the overall system model. 
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OVERALL MODEL OUTPUTS 


OVERALL EQUATION 


■ Cx 


STANDARD MODEL OUTPUT EQUATION 


^A 


C A 1 0 1 0 i 0 


—A 

*S 


1 

1 

|o 

1 

1 

|o 

1 

1 

|o 

1 

t— 

1 

0 

1 


^S 

-^OWPH 

= 

0 To 1 c ow | o 


■^OWPH 

■^OWR 


o T o I o 1 c 0H 


^OWR 

^A-S 


c A . s 1 - o,#- 



LOW ORDER SHIP MODEL OUTPUT EQUATION 


*A 


*3 

“d 


"'A-S 


—A 


-S 


TABLE 2. 3. 2.1 
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2.4 Aircraft Model 

2.4.1 Introduction 

The aircraft chosen for this study is a representative small VTOL 
vehicle with both conventional flight and vertical takeoff and landing 
capabilities. The model is based on reference ( 8 ). Figure 2. 4. 1.1 
illustrates this vehicle type. A lift fan is set in the nose with fore 
and aft, as well as lateral, thrust deflection louvers. This fan is shut 
down and doors are closed over it during conventional flight. Two more 
fans are located at the wing roots, one on each side of the fuselage. 

Each of these has a nozzle at the outlet for fore and aft thrust deflec- 

« 

tion and louvers for lateral deflection. In conventional flight,, the- 
nozzles are rotated up to produce horizontal forward aircraft motion. 

The three fans are driven by three turbine engines through shaft inter- 
connections permitting interchange of power among the fans. During hover 
(the only flight regime discussed here), all the fans are set to generate 
direct lift forces. The other forces and moments acting on the aircraft 
are due to aerodynamic and ram drag effects. Each of these effects can 
be described by nonlinear equations, and linearized for control design 
and analysis purposes. These nonlinear equations and their linearized 
counterparts are presented in the next section. Following that, the 
linearization procedure is outlined. 

2.4.2 Model Equations 

This section summarizes the original fundamental dynamic and kine- 
matic relationships, the general linearized form, and then the complete 




VTOL Vehicle Type: Fan and Louver Locations 

Figure 2. 4. 1.1 
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linearized aircraft state space model. Section (2.4.3) describes the 
details of the linearization. 

The fundamental dynamic and kinematic relationships are stated in 
» four equations: 

Newton's Law: 

Translational kinematics: 

Angular momentum: 

Rotational kinematics: 


mv_ 
— E 


^E 

4 


f„ 




at, 


T(0, V 


(2. 4. 2.1) 

'(2. 4. 2. 2) 

( 2 . 4 . 2 . 3 ) 

( 2 . 4 . 2 . 4 ) 


The force and moment on the aircraft, and m^, are composed of 
aerodynamic, ram drag, gravitational, and fan forces and moments-- The 
fan forces are commanded by the aircraft landing controller. They are 
assumed to be generated in a linear manner about their nominal values 
from the controller commands. 

The fundamental equations are linearized in Section (2.4.3). The 
results are stated below. Note that the linearization is about stationary 
hover: v_ = 0, r_ =» 0, w = 0, 9 =» 0, c * 0, and u = 0. Also note that 

— t — Ti — — D — — ~ ’ — — 

because aerodynamic forces are proportional to the square of airspeed, they 
are not small during hover over the ship deck, but they have been neglected 
to date. 

The linearized translational motion equations are stated: 

Newton's Law: 


- i t cjj Sl B + xjj£ J 


-E m 

Force summation: 

of = of + <5f 
-3 -F -R 


(2. 4. 2. 5) 


(2. 4. 2. 6) 
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Ram drag forces: 

+ ^avj, 

Fan control forces: 

%- 

Translational kinematics: 



Coordinate transformation: 

% - C E 5 % 

The rotational motion equations follow: 
Angular momentum: 

s k ‘ ^ s 2 b + h % 

Torque summation: 

6 ^= 6 ^ + 6 ^ 

Ram drag torque: 

£m_ “ T^ <Sw + cv 
^R co — B — B 

Fan control torque: 

«% - ^ M 

Rotational kinematics: 

50 = T^ 6u_ 

— to —3 


( 2 . 4 . 2 . 7 ) 


( 2 . 4 . 2 . 8 ) 


( 2 . 4 . 2 . 9 ) 


(2.4.2.10) 


(2.4.2.11) 


(2.4.2.12) 


(2.4.2.13) 


(2.4.2.14) 


(2.4.2.15) 
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The fan control force components given by Af are assumed to be 
generated as follows from the command ;u: 

Force components: 

A£ » B r Au » 5 f b t £ (2.4.2.16) 

Control state dynamics: 

c = Ac + u (2.4.2.17) 

— u— — 

The complete linearized open loop matrix differential state equation 
for the hovering aircraft can be written using the above equations. 

Table (2. 4. 2.1) presents this equation, which is of the form: 

\ = A A% + 3 A- (2.4.2.18) 

The numerical values of A and B are given in Table (2. 4. 2. 2). A math- 

ri A 

ematical block diagram of the linearized VTOL aircraft model is shown 

A 

in Figure (2. 4. 2.1). This figure presents the model structure graphically 

2.4.3 Linearization of the Equations of Motion 

The aircraft equations of motion can be divided into two parts: 
translational motions and rotational motions. In this section, these 
motions are derived as functions of fan and ram drag forces and moments. 

To date, aerodynamic effects are neglected. Next, these forces and mo- 
ments are derived in terms of fan forces, aircraft velocity, and aircraft 
angular velocity. Finally, the generation of fan forces from the controls 


is described. 
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Ram Drag Forces 
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Ram Drag Forces 
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Fan Forces 


I T M C B 
HVE 

Ram Drag Moments 


h + hA 

Internal Momentum 
+ Ram Drag Moments 


VfVt 

Fan Moments 


Euler ' Angles from 
Body Rates 


Control 

Dynamics 


Table 2.4. 2.1 
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LINEAR AIRCRAFT MODEL B MATRIX 
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0 
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1 


TABLE (2. 4. 2. 2) (Continued) 
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A general discussion of equations of motion for an arbitrary sys- 
tem is found in reference ( 9 ) . The forms and values of the matrices 
used in the aircraft model linear state equation are presented in Appen- 
dix (2. 4. 3.1). 

Translational motions are caused by forces through Newton's Law of 
motion, which is a statement of momentum conservation: 

f- = mv (2. 4. 3.1) 

'll* 

In this equation, m is the mass to which the force f„ gives acceleration 

tl 

v in the earth frame E (which is considered as inertial in this deriva- 

tion) . Linearizing this aquation and substituting in the nominal values 

for v = 0 gives: 

& 

6i E - (2. 4.3.2) 

Rewriting this equation with the force in the body frame yields the 
following equation: 

- 5 < C B + T 3 > (2. 4.3.3) 

£ 

where is the transfomration from body to earth coordinates for the 
o 

/ F 

nominal aircraft trim orientation, and T^ is the transformation between 

out of trim aircraft orientation £ and the resulting earth frame forces. 

T 

Note o9 = ( -p, 3, o) . With aerodynamics neglected, c_f p is composed of 
fan and ram drag components. df_ is given by: 

a 

3f = 6 fp + o fp ( 2 . 4 . 3 . 4 ) 
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The ram drag term <5f^ is a function of linear and angular aircraft 
velocities as derived in Appendix (2. 4. 3. 2): 

' - «r2b + < 2 

Linearized, this equation: is written 


T " 5 !2b + T V 5 ^ 


( 2 . 4 . 3 . 6 ) 


where 


T = - ^R anc * = ^cm ^ an ^ orces term given by: 



( 2 . 4 . 3 . 7 ) 


where Af is composed of the fan control force components. The kinematic 
relationship giving aircraft earth frame position r^ is very simple. The 
nominal value of r^ is zero, so the linearized equation is given directly: 

dig = dVg (2. 4. 3. 8) 

The final translational motion equation' required is a velocity transforma- 
tion: 

^B =C E^E (2. 4. 3. 9) 

Linearized about v_ = 0, this gives 

5vg = c g dyg (2.4.3.10) 

Rotational motions are caused by torques through the following 
equation:: 


h “ % 


(2.4.3.11) 
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where m^ is the torque acting on the aircraft, while Hg is the aircraft 
angular momentum. This equation is a statement of angular momentum 
conservation. Writing it in the body frame gives: 

Sg “ S B + Wg x Hg (2.4.3.12) 

The VTOL aircraft angular momentum has two parts: is the momentum of 

the engine cores and fans in the body reference frame. It is assumed to 
be constant. I^o)^ is the angular momentum of the aircraft viewed as an 

D 

entirely rigid body with inertia tensor I_ and angular velocity co_ . The 

D D 

following equation sums these two components: 

Hg - IjWg + BLj. (2.4.3.13) 

Substituting this into equation (2.4.3.12) and then linearizing gives: 

5mg = I Scog + <S<± jg x H (2.4.3.14) 

where the zero nominal angular velocity =» 0 has been inserted. 

D — 

Multiplying this equation by I_^ and introducing the notation I.. = 

o rl o 

and I. where 
A 

I A 6a)g - - Ig 1 6Wg X Hj. (2.4.3.15) 

gives the linearized result: 


% ' x m % + h 


The linearized moments o m^ are composed of fan moments Sm^, and ram drag 
moments 5m^ in the body frame. Again, aerodynamic effects, if modelled, 
would also contribute moments. The linearized moments sum to give 6 ^ : 


56 


5mg 


5 ®F + 


(2.4.3.17) 


The fan moments are a function of the control force components Af : 




m 


(2.4.3.18) 


The ram drag moments, like the ram drag forces, are generated by angular 
and linear velocities. They are given by the equation derived in 
appendix (2. 4. 3. 2): 

Sr " ‘ (2.4.3.19) 

M • • M 

Linearizing this and setting T^ = - 1^.; and T^ 11 x ) gives the 

following equation for the ram drag moments: 


‘ + <«vj, 


(2.4.3.20) 


To describe the aircraft orientation, Euler angles are used. The nominal 

aircraft orientation is heading ip = 0°, pitch angle 8 =0°, and roll 

angle <p 3 0°. Euler angle rates are integrated to determine the angles 

T ' 

which are written in vector form as 9_ = (ip, 8, <j>) . Note <56. = _8 
since the nominal 0. is zero. Euler angle rates are given by reference 
( 8 ), 

1 = T(0, ujg) (2.4.3.21) 

Linearizing this equation yields: 


U) 


" T 


<50. 


(2.4.3.22) 
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This completes the description of the aircraft motion model given 
the fan forces. The description of the command and actuator, coordina- 
tion and dynamics follows . 

The generation of fan forces and moments from commands is accom- 
plished in four steps. The command _u is passed through simple dynamics 
to give the command state vector _c. This is transformed into actuator 
variables represented by Au, which in turn are transformed into the force 
components of the fans, Af . Finally the total fan forces, f^., and 
moments, m^, are obtained by a pair of transformations acting on Af. An 
expansion of this description follows. 

The command u has six components representing the co mm anded inputs 

to the x,yv2,9,<j), and motions. In component form it is written as 

T 

_U = (5x, 5y, 5z, 5<j), 59, Sip) . The dynamics given the command states 

are uncoupled and first order. For example, 5x c is the state of a first 
order dynamical system whose only input is 5x. The overall command state 
generator system is described by the equation: 

c - A c + u (2.4.2.23) 

— u — — 

T 

where h = (<$ x c » 5y c> ^ z c » ^^c^ ' * These dynamics are used 

to model the accumulated dynamics of the control system, control actua- 
tors, and thrust generating mechanism. A u is chosen as a zero matrix so 
the open loop command states are the integrals of the command components. 
The next step in the generation of the forces and moments is the trans- 
formation from command states to actuator variables. This transformation, 
shown in the next equation, is chosen to yield uncoupled aircraft motions 
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from the uncoupled motion, commands : 

Ay = Be. (2,4.3,24) 

Au contains longitudinal and lateral thrust deflection angles, and thrust 

magnitude deviations for each of the three fans. The quantities are 

T 

Au *» (Ac^, Aa 2> Aa 3 , A0^, A0 2 , A0 3 , AF^, AF G2 » AF G3 ) • b t transforms 
command states into actuator variables by the following scheme. A for- 
ward motion command state is converted to equal aftward deflections 
in each lift/cruise fan nozzle and in the forward lift fan louvers. 
Likewise, lateral and vertical command states <5y and <5z are converted 
to lateral thrust deflection and thrust increase respectively. Roll 
is generated by equal and opposite thrust level variations in the lift / 
cruise fans. <5Q c> the pitch command state, is transformed into increased 
lift fan thrust and decreased lift/ cruise fan thrusts in amounts causing 
rotation and not translation. Yaw command state 6 i^ c is similarly trans- 
formed into yawing moment by deflecting the lift fan thrust to one side 
and the lift/cruise fan thrusts to the opposite side. The following 
equation describes how the actuator variables Au transform into fan force 
components: 

Af - b f Au (2.4.3.25) 

where Af = (AF xl> AF^, AF x3 , AF yl> AF y2 , AF y3> AF^, AF^, AF^) 

gives the force components for each of the three fans. These fan com- 
ponents are combined to form the total forces on the aircraft by: 
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5f F = ijM (2.4.3.26) 

T 

where 6f^ = (X^, Y^,, Z^) . Similarly, the moments are given by: 

6m F * T^Af (2.4.3.27) 

T 

where iSm^, = (L^, M^, N^) has roll, pitch, and yaw moments respectively 

as elements . 

2 .5 Wave Model 
2.5.1 Introduction 

The purpose of the wave model is to generate the stochastic vari- 
ation of sea surface height at some point along the ship’s length. This 
is used as the disturbance driving the ship motion model. The use of 
sea surface height to drive ship motion is a commonly used approximation 
(10. ) . It is a major simplification of the distributed hydrodynamic 
forces acting on a ship hull, but it does result in an overall wave-ship 
model whose response is characteristic of actual ship motion. General 
wave surface and ship motion modelling are described in reference (10 ) . 

The Neumann ideal power spectrum for wind generated waves is chosen 
as the starting point for the wave modelling. Reference ( 11 ) gives 
this spectrum. However, for control design and systems analysis, a 
linear time invariant dynamic system driven by a gaussian white noise 
and written in state space form is desired for the wave model. Therefore, 
the result of the wave modelling is such a state space representation 
having an output spectrum approximating the Neumann spectrum. 
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The sea condition, and ship speed and heading are required to 
specify the wave spectrum encountered by the ship. Ship information is 
needed in order to transform the wave spectrum from that seen at a 
stationary point to that seen from the moving ship. The sea and ship are 
given specifications which provide a realistic, challenging situation 
for the aircraft landing controller to operate in. The specifications 
are for a "sea state five" (defined in reference (12 )), a ship heading 
such that the relative wind over the deck is undisturbed by the ship's 
superstructure, and a ship speed compromising low wind speed over the 
ship deck and high through the water ship speed. The wind speed gene- 
rating the sea state five ocean waves is set at 16 knots ( IQ ) . In 
steady state over a wide stretch of ocean where the waves will not be 
fetch limited, the wave heights will reach over 6.5 feet about 30% of 
the time. (This is generally called the significant wave height, 
measured from trough to crest.) The corresponding period of these waves 
will be 7 seconds (10). Ship velocity through the water is chosen so 
that the wind sweeps the ship landing pad at not more than 20 knots, 
and at a heading (9^) near 45°. This angle avoids air turbulence over 
the pad due to the ship superstructure. The specification of a low 
relative wind speed keeps aircraft aerodynamic forces small; and, hence, 
helps to justify neglecting those forces. In addition, the ship should 
be cruising at a maneuverable speed. For these reasons, the ship 
velocity is chosen at about 10 knots, and at an angle (0^) of 70° into 
the wind and waves . 

The following section describes the ideal wave spectrum and the 
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final model. The derivation of the latter is outlined in the subsequent 
sections . 


2.5.2 Model Equations 

This section gives the form of the ideal Neumann spectrum for both 
fixed and moving observers. Then the state space approximation adopted 
as the stochastic wave model is presented. The following two sections 
describe how this model is obtained from the ideal spectrum. 

The Neumann power spectrum for wind generated ocean waves has the 
following general form as a function of wave frequency 


SCu^) = a 2 0 J w s exp {-Ibu^ 2 } 


( 2 . 5 . 2 . 1 ) 


where a and be are constants depending on the wave height and period 
characteristic of sea conditions as explained in the next section. Trans- 
forming this spectrum into the one observed from onboard the ship moving 


with speed v_ and heading 9„ into the waves gives S (oj ) : 

S w e e 

S e (w e ) = a 2 ui^ 6 exp { -2ba^ 2 }[l - ( 4v s w e cos9 w /g ] 


where: 


1/2 


- 1/2 


( 2 . 5 . 2 . 2 ) 


-1 


= { -1 + [ 1 + (4v s u e cos6 w )/g] } [(2v g cosO w )/g] 

(2.5.2. 3) 


Note g is the acceleration due to gravity. In figure (2. 5. 2.1), the 

1/2 

amplitude spectrum (S^Cas^)) is plotted for comparison with the finite 
order state space model. Note that due to the exponential factor the 
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ideal spectrum drops off faster than any finite order model spectrum at 
low frequencies. 

The finite order model was chosen so that its amplitude spectrum 

1/2 

would match the ideal amplitude spectrum, (S (w )) , closely. This 

model was formulated with the constraints of finite order, simple struc- 
ture, and gaussian white noise input. The result is shown in Table 
(2.5. 2.1). The amplitude spectrum of this model is plotted in Figure 
(2. 5. 2.1). For the decade of frequencies centered at the peak amplitude, 
the ideal and model spectra match well. Outside this range the spectra 
diverge, but the divergence is only significant at low and high fre- 
quencies where the amplitudes are very small. 

The output of the model, y , is extended in the ship modelling 
sections, where it is used as the ship model input. Ship pitch and 
heave inputs are related by phase, which can be accounted for in the 
output vector of the wave model; however, no phase relationship of roll 
with pitch or heave was determined, so the roll input is generated by a 
second, identical, independently driven wave model. 

In the next two sections, the state space model development is 
described in more detail beginning in the following section with a more 
detailed explanation of the Neumann ideal wave spectrum. 

2.5.3 Ideal Ocean Wave Spectrum 

Various spectral densities for ocean waves are available to naval 
architects in frequency domain analysis of ship motions (10, 11 ). Of 
these spectra, the theoretical Neumann spectrum is chosen here. It is 
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OCEAN WAVE STATE SPACE MODEL 

^OW = A ow^ow + B OW P OW 
7 0W = G 0#0W 



TABLE (2.5. 2.1) 


Continued 
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OCEAN WAVE STATE SPACE MODEL (CONTINUED) 


X 1 


0 

1 

0 

0 0 

0 


X 1 


0 

X 2 


-.50126 

-.92040 

0 

1 0 

0 


x 2 


0 

X 3 


0 

0 

0 

1 0 

0 


X 3 

U. 

0 

X 4 


0 

0 -. 

50410 

-.92300 0 

1 


x 4 


0 

X 5 


0 

0 

0 

0 0 

1 


X 5 


0 

X 6 


0 

0 

0 

0 -.49844 

-.91780 


x 6 


.12810 

• 


y - C -8.0555 -14.750 0 15.980 0 0 ] x„ rT 

— ow 


TABLE (2. 5. 2.1) (Continued) 



66 


simple, yet approximates actual sample power spectral densities reasonably 
well. As presented in reference (11), it has two parameters, and 
which characterize average apparent wave period and minimum height of the 
1/3 largest waves, respectively. (In the modelling, is assumed to be 
the period of waves of height The Neumann expression for the wave 

power per unit frequency, as seen by a stationary observer is: 

s( “w’ ■ "5/3 I » (3,!,|3,I W 2i) " 6 

(2. 5. 3.1) 

Travelling on the ship through the waves, a transformed or distorted 
spectrum is observed. The transformed spectrum can be obtained using a 
Jacobian transformation. This transformation comes about in the fol- 
lowing way. A moving observer sees the same physical waves, as a sta- 
tionary observer, but at a different frequency and power density, say 
S e (io e ). The new observed frequency is called the frequency of encounter. 
Since both observers must see the same total wave power, the follox^ing 
equality holds: 

Q f S(w w )dw w = Q f S e (w e )du) e (2. 5. 3. 2) 

The Jacobian transformation, or change of variables formula (13 ), gives: 

S (u ) - S(ui ) (2. 5. 3. 3) 

6 6 • 3u> 

U e “ f( V 
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To solve this expression for S e (w g ) , the frequency tranformation, 

U = f (ojy) , given in reference (10 ) is required. This reference 
derives f(o^) from the properties of deep water waves and the frequency- 
wavelength-velocity of observer doppler shift relationship. The result 
is: 

u e = f (a^) - oj w cos 0 w ( • (Vg/g)^ + l/cos9 w ) (2. 5. 3. 4) 


or 


-1 + (1 + Ao^VgCosQ^g) 1 ^ 2 

W 2v g cos0 w /g 


(2.5. 3.5) 


where oo^. is the frequency observed by a stationary observer, v^ is the 
ship speed, 0^ is the ship heading with respect to the direction of the 
wind and wave source, and g is the acceleration due to gravity. The final 
expression for is: 


,2 ok 4v cosq.. -1/2 

Se'Cwe) = ^ ex P [1 + — — ~ “e 1 (2.5. 3. 6) 

where 

a = H 1/3 (3/128 tt s ) 3/4 T~ 5/2 (2. 5. 3. 7) 

b => 3 it 2 T^ 2 (2. 5. 3. 8) 

and is given by equation (2. 5. 3. 5). This is the power spectral den- 
sity used to determine the linear finite dimensional stochastic wave 
model in the next section. 

An alternative derivation of S e (tu e ) can be accomplished based on 
the Taylor "frozen field" hypothesis used for wind gust modelling (14 ) . 
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The frozen field hypothesis assumes that the spatially' ’distributed gust 
process W(x,t) is a fixed field in time which moves spatially at velocity 
v . Hence, it has the property: 

W(x,t) = W(0, t— x/ Vp) (2. 5. 3. 9) 

For ocean waves, the hypothesis would be applied to individual wave 
frequency components, which may be thought of as fixed in time, but 
moving spatially with frequency-dependent velocity v^Cw) . 

The next section explains how a state space model of the ocean waves 
is derived, whose spectrum approximates the Neumann spectrum. 

2.5.4 Model Derivation 

The purpose of this section is to describe the derivation of a state 
space ocean wave model having approximately the same spectrum as the 
Neumann spectrum in the previous section. 

The Neumann spectrum of equation (2. 5. 3. 6) represents an infinite 
order linear system. The reason is that it contains an exponential 
factor, which at low frequencies causes the spectrum to drop off faster 
than that of any finite order system. On the other hand, it is desirable 
to have a simple, low order model for computations and analysis. The 
construction of such a model follows. 

Modelling is a system identification problem: a structure must be 

hypothesized, and its parameters estimated. The hypothesized structure 
is a gaussian white noise driven single input-single output (SISO) system. 
Its dynamics consist of three nearly identical second order damped 
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oscillators in cascade. This structure has the advantage that it con- 
tains only three parameters affecting the output, while satisfactorily 
matching the ideal spectrum of the previous section. The small variation 
in the parameters among the three second order systems is to aid numerical 
algorithms. The model transfer function is: 


Ocean surface height 
noise input 


(s) 


Cs 


1 4 


n (s 2 + 2£u) .3 + w 2 . ) 
i=i 01 


(2. 5. 4.1) 


where s is the Laplace variable, C is the root locus gain, u!> ^ is the 
ith second order system natural frequency, and £ is the common damping 
ratio of the cascaded second order systems. The state space differential 
equation for this is written in Table (2. 5. 2.1) . 

Having chosen a hypothesized structure, its parameters need to be 
estimated. This is accomplished by matching the square roots of the 
power spectral densities of the finite order model with the infinite 
order ideal reference spectrum. If y is the output of a single output 
system with power spectrum S(u)), then the quantity under consideration is: 

(S(u)) 1/2 = C&E [ (y (t) - m(t))(y(t+T)-m(t+T)]} (<o)) 1/2 (2. 5. 4. 2) 

where 

m(t) = E[y(t) ] (2. 5. 4. 3) 

and (w) is the Fourier transform, while E is the expectation 

operator. For single input-single output linear time invariant systems 
of the model form given in Table (2. 5. 2.1), the square root of the power 
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spectrum is: 

" I C OW^ uI “ A ow ) B ow^ ^ sp ow (co)) 

( 2 . 5 . 4 . 4 ) 

where Sp^^Cw) is the power spectrum df the input P 0W^^ * Since Pgy(t) is 

1/2 

white noise, Sp^^Cw) is a constant, independent of frequency. (S (ui) ) 
is refered to as the amplitude spectral density. 

Four specific features of the ideal and model amplitude spectra are 
matched. These are: peak amplitude and frequency, peak width, and roll 

off in the first decade on each side of the peak. The roll off towards 
low frequencies is obtained by four zeroes at the origin of the root 
locus plot. Four poles located at the ideai amplitude spectrum peak 
cancel the effect of the four zeroes at high frequencies, while two more 

poles at the same location produce the desired high frequency roll off. 

/ 

Thus, the attenuation away from the peak is determined by the model 
structure itself. Since the peak frequency for a second order system 
with a zero at the origin is near its natural frequency 0)^, the values of 
w q 2 > u n 3 are chosen near the ideal spectrum peak at 0.708 radians/ 
second. The values used are: 0.706, 0.708, and 0.710. The damping 
ratio is set so the width of the model and ideal amplitude spectra 
roughly match over the first decade of drop from the peak. The value 
used is £ = 0.65. For comparison, the ideal Neumann and the model 
amplitude spectra are presented in figure (2. 5. 2.1). Verification of 
the model validity is by the close match illustrated in this figure. A 
further check is made in Section (2.6.5) where a simulation ! is presented 
to verify the overall wave-ship model behavior. 
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The wave model parameter fit was made by eye judgement; however, a 
better fit might be obtained with a more precise method such as minimum 
least square error fit. In view of the fact that the ideal wave spectra 
are themselves only rough approximations of wave behavior, more elaborate 
fitting methods were not deemed to be justified. 

2.6 Ship Motion Model 
2.6.1 Introduction 

This model represents the dynamics of the ship by transforming ocean 
surface height into ship motion response. In particular, the landing 
deck motions are represented. The model has inputs for the roll and 
pitch-heave dynamics of the ship, and outputs for the resulting orienta- 
tional and translational landing pad motions. Fore-and-aft sway, yaw and 
lateral sway are not modelled since it is assumed that their responses 
are smaller than the others, and so less Important to the aircraft 
landing problem. The ship modelled is a representative small ship like 
those of the Allen M. Sumner Class. Ships of this class are described 
in reference ( 15 ) . . They are approximately 400 feet in length and dis- 
place roughly 3000 tons. As described in the wave model introduction, 
the ship is assumed to be cruising at -10 knots, with a heading angle of 
70° into the waves. The waves are generated by a 16 knot wind. Ship 
motion magnitude parameters are sized to produce typical ship roll, 
pitch, and heave amplitudes of 5 degrees, 2 degrees, aiid 5 feet, res- 
pectively, as suggested by reference ( 7 ) . The heave amplitude cor- 
responds to about one foot of heave per foot rise in the sea surface; 
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however, roll and pitch do not appear to correspond with a similar static 
relationship using wave slope. 

The next section describes the linear, time invariant ship model and 
its origins. The other sections describe the identification of dynamics, 
the calculation of outputs, and the verification of the overall wave-ship 
model. 

2.6.2 Model Equations 

i 

This section presents a brief review of the ship motion modelling, 

9 

and then gives the model results. 

The ship motion analysis techniques used in the modelling are found 
in reference ( 10 ) and in reference ( 16 ) . The second is the classic 
paper in the field. In these references, the basic assumption is that 
the overall ship response is the superposition of responses to the fre- 
quency components of the ocean waves. The important result is that to 
a good approximation, the ship response to irregular seas can be repre- 
sented by a cascade of the wave model and a shaping filter, where the 
filter represents the dynamic response of the ship to the waves. Another 
result is that longitudinal and lateral ship dynamics may be treated as 
uncoupled for most theoretical and experimental analyses . For the ship 
modelled here, coupled pitch and heave dynamics data is available from 
a towing tank experiment ( 17 ) , However, data on rolling motions has 
not been located. This lack of roll data is overcome by observing and 
modifying typical rolling data for other ships ( 18 , 19 ). Since roll 
dynamics are not -used in the control system design and analysis, the 
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purpose of including roll is to generate a complete model of the major 
ship responses. Recall from the introduction that surging, yawing, and 
swaying are not expected to be as important as pitch, heave, and roll 
responses, and so are not modelled. 

Characteristics of the shaping filter representing ship response 
are used for the identification of the state space ship model parameters. 
In this linear time invariant model, pitch, heave, and roll modes are 
each represented by a separate second order approximation of a pass band 
filter. ‘To generate the desired phase between pitch and heave, the heave 
input is lagged behind the pitch input. Since the ship inputs are wave 
model outputs, the single wave model output is augmented to include a 
second phase lagged quantity. The two wave output quantifies used for 
ship inputs are y pw and y^. Since no roll-pitch phase relationship 
could be determined, roll is driven by a separate identical, but statis- 
tically independent wave model. The lack of roll information can be 
attributed to two factors. Rolling motions are not as critical as 
pitching and heaving in determining a ship’s ability to maintain heading 
and speed in rough seas; and in those cases where rolling is unacceptable, 
it can often be reduced to acceptable levels, with the use of roll damping 
devices . 

The end result of the modelling is shown in the tables . The ship 
model equation with the three independent second order pass band filters 
representing roll R, pitch P, and heave H dynamics is shown in Table 
(2. 6. 2,1). In this equation, is the natural frequency and £ the damp- 
ing ratio of the indicated dynamics. y R ^, y pw , y^ are the three ship 
inputs from the wave models. Next, the output equation is given in Table 
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SHIP STATE SPACE MODEL DYNAMICS 
(The output equation follows in Table 2. 6. 2. 2)) 

x s - A s Xg + B s Iwsc 
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TABLE (2. 6. 2.1) 


(Table continued on next 
page) 
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(continuation of Table (2. 6. 2.1)) 


*R1 


*H2 


*P1 

__ 

*P2 


*H1 


*H2 



0 

0 

0 

0 

0 

0 

-0.10693 

-0.02616 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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0 0 0 
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0 0 0 
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0 0 0 
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TABLE (2. 6. 2.1) continued 
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(2. 6. 2. 2). (x^p, y ^p , z ^-p ) is the position of the landing pad in 
earth coordinates, while ( <P S , 0g., are the Euler angles of the ship 
orientation. In this equation, C is generated by the matrices 
Sgi» S*^, Sg^, and Sj^. These are described in Section (2.6.4). The 
coupling structure between the wave output and ship input is presented 
in Table (2. 6. 2. 3). x„ r _„ and x_ r _ are the state vectors of the wave 
models associated with the ship pitch-heave and roll dynamics inputs, 
respectively, w ^ an ^ ?£ are a ^ so f rom the wave model. Finally, the 
overall wave-ship model is written in Table (2. 6. 2. 4). In the overall 
structure, note that P 0W p H and PqwR are in ^ e P en ^ ent gaussian white 
noise sources driving the two wave models. 

In the next section, the ship dynamics identification is described. 


2.6.3 Ship Dynamic Modelling 

This section discusses the identification of the structure. and 
parameter values for the ship model. The discussion begins with the 
important analysis references and sources of data. Then it proceeds to 
the model structure and parameter estimation. 

The modelling approach for the ship was to take the frequency 
domain response representation of references ( 16 ) and ( 10 ) and approx- 
imate it in the time domain by a low order linear system. The ship 
motion data for estimating the model parameters came from several sources. 
For the ship modelled, as for many ships, only pitch-heave information 
was available. For roll, information on other ships was employed. Data 
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SHIP MODEL OUTPUT EQUATION 


*s = C s^s 



TABLE 2. 6. 2. 2 
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WAVE OUTPUT - SHIP INPUT COUPLING 
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C QW = [ -Go*, -2 Go 2 C 2 , 0, C, 0, 0 ] 

C PH * C -qC<0 2 , -2qOo 2 C 2 + qCb, 0, qC, 0, 0 ] 

C QW = [ -8.0555, -1.4750, 0, 15.980, 0, 0 ] 

Cp R = [ -4.0696, 3.2535, 0, 8.0731, 0, 0 ] 


TABLE (2.6.2. 3) 
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OVERALL WAVE-SHIP MODEL STRUCTURE 
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TABLE (2. 6. 2. 4) 
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from the following sources was 'used for the parameter estimation. Experi- 
mental towing tank pitching and heaving data for a model of the 383 foot 
long Allen M. Sumner Class destroyer was obtained from reference ( 17 ) . 
Data was also available from model tests of a similar 500 foot ship in 
reference ( 20 ) . Full scale rolling data for a 364 foot destroyer was 
found in reference ( 19 ) . In reference ( 18 ) , model data was found for 
the pitching, heaving and rolling of an aircraft carrier. 

The model structure is;. now familiar. It is a compromise between 
simplicity and realism, while taking into account the available data. 

The response of the ship to ocean surface height is assumed to be linear. 
Longitudinal and lateral dynamics are decoupled; and, in fact, pitch and 
heave are decoupled, except for a relative phase relation at their inputs. 
Qualitatively, the response of the ship to waves is roughly like a pass 
band over a narrow frequency range. In fact, the three modes behave much 
like highly tuned second order systems. The references (10), (17) and 
(21) illustrate this characteristic. In these references, the response 
amplitude operator is plotted. This is the square of the magnitude of a 
system's amplitude response. Based on this observation, each mode of 
response is hypothesized as a second order system. The overall ship model 
is thus three independent second order systems driven by individual in- 
puts as shown in the previous section. Table (2. 6. 2.1). The ship state 

i 

vector, Xg, is given: 

Xg = (x^» Xp 2 , x H1 , x 02 ) (2. 6. 3.1) 

where for each mode R, P, and H, the first state is proportional to the 
orientation or height of the ship in the earth frame, and the second 
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state is the derivative of the first. Magnitude scaling is accomplished 
at the output. 

Pitch and heave parameters were estimated first. The parameters are 
the natural frequency and damping ratio of each mode. Reference ( 17 ) 
was the primary source of data here. The first step was to plot amplitude 
response as a function of frequency of encounter. For a given speed, with 
the ship heading directly into regular waves, the motion amplitude was 
recorded for various wavelengths. This information from the report can 
be plotted as motion amplitude versus wave frequency. A transformation 
yields a plot of amplitude versus frequency of encounter for the given 
ship velocity. The transformation from to to is the same one used in 
the wave model derivation. 

For the pitching case, it is easy to find the resonant frequency to 
use as the natural frequency. A value of a> n p =» 0.65 radians /second was 
selected. By comparing the attenuation of the response magnitude at two 
frequencies, a damping ratio can be determined. The value = 0.25 
was selected. 

For heave, the data does not show a peak. From reference ( 20 ) , 

however, information was available for a 500 foot destroyer. There, the 

heave and pitch resonant frequencies were the same. On this basis, go 

tin. 

was chosen equal to cu^-p . The heave mode damping ratio was set close to 
the pitch damping ratio since the amplitude responses have about the 
same attentuation rate near the resonant frequency. The model value is 
Cg - 0.24. 
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A pitch -heave phase relationship was given in reference ( 17 ) . This 
was used to generate separate pitch and heave mode inputs to the ship 
from the wave model output. As shown in Table (2. 6. 2. 3), this was accom- 
plished by giving one wave model two outputs, y ?w - and y^. (Note the 
second wave model drives the roll mode.) The phase relation is as follows 
For very low frequency waves, that is ones with long wavelengths, pitch 
will lead heave by 90 degrees of phase. This is because pitch will be 
at the maximum positive value when the ship center of rotation heaves 
upwards through the mean position as the ship climbs a wave to its peak. 

l 

At the crest of the heave motion, the pitch attitude will return to zero, 
and so on. At a wavelength corresponding to the frequency of encounter 
equal to about 0.68 radians /second, the data shows that pitch leads heave 
by only approximately 60 degrees. These two phases . suggest the following 
frequency domain relation between 7p^» the wave output driving ship pitch, 
and y^, the output driving, ship heave: 

y HW = <1 ( s + b ) S_1 7 PW (2. 6. 3. 2) 

where q - 0.51 is set to give the heave and pitch inputs the same magni- 
tude at pitch and heave natural frequency, and b = 1.3 gives the desired 
phase relation. The equation for y^ in the time domain is given in 
Table (2. 6.2.3) . 

Roll motion data for the Allen M. Sumfter Class was not found; how- 
ever, data from other ships was used as a guide for modelling a second 
order system to it. From reference (18) and other reports, it was 
noted that the rolling frequency of a ship is about half the pitching 
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and heaving frequency. Thus, = 0,33 was selected. From the same 
report, an est ima te of the roll damping ratio was made. The ship was 
an aircraft carrier with a length to breadth ratio of 7.7, and its 
damping ratios were ■ 0.26, fp ** 0.20, and = 0.032. These values 
roughly fit the ratio: 


length ^ ^P 


beam 


’R 


( 2 . 6 . 3 . 3 ) 


Assuming this is also approximately true for the Sumner Class destroyer, 
where the length/beam ratio is 9.4, the value 0.027 would be chosen for 

Another report, reference ( 19 ) > for a 364 foot destroyer gave data 
from which its damping ratio was calculated at £_ = 0.05. As a compro- 
mise between these two values, the damping ratio for the model was set 
at = 0.04. 

Nb. « phase relation: 1 is known for roll with respect to pitch, so the 
roll mode of the model is driven by an independently excited copy of the 
wave model, as stated before. The wave model output exciting roll is 
designated y^ as shown in Table (2. 6. 2. 3). 

In the next section, the output equations and the scaling of the 
motions are discussed. 


2.6.4 Ship Model Outputs 

This section describes the outputs of the ship model. 

The ship responses desired in the output are those of the landing 
pad (LP) in earth coordinates. The output vector includes the landing 
pad position, velocity, orientation and angular rate. It is written: 
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^ X ELP* y ELP* Z ELP * * 9 S* ^S* 

®S* ^ 


X ELP* y ELP’ Z ELP* ' 

(2. 6. 4.1) 


where (^p, y ELp , 2 CTP ) is the position and (<J> g , '0g, i//g) are the Euler 
angles of the landing pad orientation. Note that ip$ =» 0 and \p g ** 0, 
since yawing motion is not modelled. 

Tile ship state vector must first be scaled to the correct magnitude 
before the output is computed. The desired magnitudes are expressed as 
root mean square (fms) values. They are 5 degrees, 2 degrees and 5 feet 
respectively for roll, pitch and heave. To achieve these values, the 
scale factors AMPRLL, AMPPCH, and AMPHVE are calculated to give the cor- 
rect corresponding diagonal elements of the steady state,, state covariance 
matrix. The values are shown in Table (2. 6. 4.1). There, the matrix 
scales the roll, pitch and heave states, while S g v. scales the rates. 

The motion at the landing point is found by multiplying the scaled 
states by T^gg, which is given in Table (2. 6. 4. 2). To retain the pro- 
perties of linearity, T £ g£ is evaluated at the nominal orientation of 
the ship, 0 g =0, <j> g - 0. The quantities x gLp , y gLp , and Zg Lp in T £SC 
are the coordinates of the landing pad in ship coordinates. The values 


are -175.0, 0.0, and -16.0 feet respectively. 


Xg £p was chosen about 


half the length of the ship aft of the center of rotation, while the 


vertical deck position z gLp is estimated at 16.0 feet above the center 
of rotation. (Note vertical [distances are positive downwards in earth 


and ship coordinate systems.) y gLp is set t0 zero » directly on the 
ship center line. 
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SHIP STATE SCALING MATRICES 

f 

Sg^ Scales roll, pitch and heave states to radians and feet. 

S^ Scales roll, pitch and heave derivative states to 

radians/second and feet/second 



AMPRLL 

0 

0 

0 

0 

0 

S S1 - 

0 

0 

AMPPCH 

0 

0 

0 


0 

0 

0 

0 

AMPHVE 

0 


0 

AMPRLL 

0 

0 

0 

0 

S S2 " 

0 

0 

0 

AMPPCH 

0 

0 


0 

0 

0 

0 

0 

AMPHVE 


0.01821 

0 

0 

0 

0 

0 

S S1 - 

0 

0 

0.01 

0 

0 

0 


0 

0 

0 

0 

1.261 

0 


0 

0.01821 

0 

0 

0 

0 

S S2:: - 

0 

0 

0 

0.01 

0 

0 


0 

0 

0 

0 

0 

1.261 


TABLE (2. 6. 4.1) 
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COORDINATE TRANSFORMATION FROM SHIP STATES TO SHIP LANDING PAD POSITION 

IN THE EARTH FRAME 

Abbreviations: S represents sine 

c represents cosine 
The subscript. S on 0 and 4> is dropped. 


s0c * y SLP + s6s * Z SLP - s6x SLP + ' c6s4,y SLP " c9c<J)Z SLP ° 

T ESC = - S0y SLP + C ^ Z SLP ° ° 

c0ct f> y SLP + c9st f >2 sLP - C0X SLP ' S0S ^SLP “ sQ c< ^ z slp - 1 

nomi- 
nal 
0=0 
<j) = 0 



0.0 

-16.0 

0.0 

T = 

ESC 

16.0 

0.0 

0.0 


0.0 

175.0 

-1.0 


TABLE (2. 6. 4. 2) 
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■<; The output matrix Cg is given in Table (2. 6. 2. 2). The first three 
outputs are the landing pad position, formed by first scaling the 
position states by Sg^, and then transforming them to give the pad 
coordinates. Similarly, the rates are obtained by scaling the rate 
states by S^, and then transforming them to give pad position velocities. 
The second and fourth triplets of outputs are the angular orientations 
and. rates. Sg^ and Sg 2 do the scaling. However, the third element for 
height is zeroed, so the third element represents yaw orientation or 
rate; hence the notation S*^ and S* 2 . 

2.6.5 Overall Wave-Ship Model Verification 

Qualitative features of the ship motion model behavior are verified 
in this section. 

Typical time histories of ship motions in seaways are available for 
various ships and sea conditions. A typical pitch motion trace for a 
ship in a seaway is taken from reference ( 16 ) , and redrawn in 
Figure (2. 6. 5.1). No time scale was given for this trace in the reference. 
Nevertheless, two characteristics may be noted; First, the period is 
very regular; and second, the pitch amplitude varies quite irregularly. 

The regularity of the period corresponds to the concentration of power 
in the ship motion power spectral density almost entirely at the natural 
frequency. For comparison, a segment of a simulation of the wave-ship 
model is also shown in the figure. Qualitatively, the simulation trace 
displays the same two characteristics as the actual ship motion. 
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Comparison of Pitch Motion of an Actual Ship in a Seaway and the Ship Model Simulation 
Figure 2. 6. 5.1 


89 


2.6.6 Low Order Ship Motion Model 


This section describes a simpler model for ship motion. In the hover 
control chapter, it is used in a comparison demonstrating the effect of 
modelling accuracy on hover control performance. Historically, this model 
preceded the separate cascaded wave and ship models of the previous sec- 
tions. Those models were formulated when the low order model was found 
to be an insufficiently accurate representation of typical ship motion. 

The problem with this model is that it provides a relatively broad pass- 
band of motion energy, while true ship motion actually occurs in a quite 
narrow passband. A comparison between the two models is shown in Figure 
( 2 . 6 . 6 . 1 ) . 

In the low order model, each mode of ship motion (designated by i) , 
is assumed to be driven directly by an independent, zero mean, gaussian, 
white noise of intensity 1.0. The response i to input i is second 
order with transfer function: 


response i , ^ 

-r “ 


input i 


(s : 


§-r s 


+ 25 iV + 


“il> 


( 2 . 6 . 6 . 1 ) 


The natural frequency is chosen as 0.52, 0.97 and 1.0 radians/second 
for roll, pitch and heave, respectively. These values are according to 
references (22) and ( 7 ). The gain g^ is set so the root mean square 
output amplitudes would be 5 degrees, 2 degrees and 5. feet, as done in 
reference ( 7 ) and in the standard model. The damping ratio is set 
to a small value, 0.05, to give a highly tuned second order response. 

The resulting algebraic form is identical to that for the regular ship 



Ship Motion Model Comparison 


Figure 2. 6.6.1 
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model shown in Table (2. 6. 2.1), with the exception that the input is no 

T T 

longer (y RW > y pw , y^) , but (p R , p p , p R ) . The output equation is also 
the same. The final low order ship model dynamics are shown in Table 
(2. 6. 6.1). Since the output equation is the same as in Table (2. 6. 2. 2), 
it is not repeated here. 
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LOW ORDER SHIP MODEL DYNAMICS 



TABLE (2. 6. 6.1) 
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CHAPTER 3 
HOVER CONTROL 


3.1 Introduction 

In this chapter, a hover controller design for the VTOL aircraft 
vertical and pitch dynamics is developed. The purpose of the controller 
is to hover the aircraft at a constant vertical distance above the 
landing pad of the ship. Both of these vehicles are described in the 
modelling chapter. The challenge to the controller is the random motion 
of the ship, which cannot be predicted exactly. For this problem, as it 
will be defined in this chapter, the derived design is optimal. There is 
no other design solution whose performance would be better. However, 
this does not exclude a redefinition of the problem to obtain another 
more desirable design. A second problem and approach are in fact recom- 
mended here, and investigated in the next chapter. The recommendation 
comes about because the design solution of this chapter is representative 
of the optimal design under idealistic assumptions, and could easily 
become unacceptable under realistic conditions 

This chapter has three main parts. They are: the design goals, 
the method to achieve the goal, and the results of the hover controller 
design. The first part explains the goals of hover control, including 
the basic assumptions and the theoretical tools available for the con- 
troller design. The second part explains how the assumptions and theory 
are employed in the design method. Finally, the third part discusses the 
results and how well they achieve the desired goal. 
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3.2 Control Problem Definition 

3.2.1 Design Goals 

The design goal is to determine an aircraft, shipboard landing 
technique based on modern control theory. An important constraint is 
that the technique should require only a small amount of thrust control 
authority. In this chapter, the goal Is posed as a tracking problem. 

The aircraft is to hover at a constant separation above the ship deck, 
as illustrated in figure (3.2. 1.1). In this task, it is desired to 
minimize the deviation of the aircraft from its desired tracking position. 
The steady state linear quadratic gaussian (SSLQG) regulator gives the 
optimal solution to this tracking problem when both the aircraft and 
ship equations of motion are linearized about their nominal operating 
point; and their separation is defined as an output for which a quadratic 
cost function is chosen. The SSQLG regulator solution provides a steady 
state control law which minimizes this cost in the presence of gaussian 
white plant noise. No other control can reduce the cost further. 

3.2.2 Assumptions 

Certain assumptions in addition to those made in the linearized 
model are useful. The model gives linear equations for the aircraft 
dynamics near stationary hover, and for the ship moving at a constant 
velocity in ocean waves. To simplify the control design and analysis, 
the following additional assumptions are made. 

The first assumption is complete state measurement. This not 
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only simplifies the control dynamics, but is an important key to the 
interpretation of results as being representative of the "best 
theoretically achievable" performance. Under practical circumstances, 
measurements are made with noisy sensors of quantities related to, but 
not necessarily identical with, the ones of interest. For feedback 
control, the elements of the system state are the desired quantities. 
Thus, an estimator such as a Kalman filter is generally required to 
estimate the quantities. The form of the estimator depends on the. 
measurements available. In general, it introduces additional dynamics 
into .the overall system, retarding the overall system response. By 

i 

assuming that uncorrupted state measurements are directly available, 
the estimator is not necessary and the best, possible performance is 
achieved. Results presented here should therefore be interpreted as 
lower performance bounds not necessarily attainable in practice. 

The second assumption concerns lateral-longitudinal decoupling. 

For the ship model, this assumption is justified in Chapter Two, where 
the ship is represented as two separate decoupled systems, with roll 
motions entirely separate from pitch and heave motions. As noted in 
Chapter Two, this is common practice in ship motion modelling. The 
aircraft model, on the other hand, is not constructed with separate long- 
itudinal and lateral directional dynamics (although that is also common 
procedure for conventional aircraft ( 9 )). Indeed, the linear aircraft 
model exhibits some coupling terms between the two types of dynamics. 
However, these terms are quite weak and can be neglected for most 

early analysis. 
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The next level of assumption is that forward horizontal aircraft 
dynamics do not couple into vertical translation and pitch orientation 
dynamics. This may also be verified by refering back to the aircraft 
model. Hence, the aircraft is left with only vertical and pitch motion 
dynamics for the design and evaluation. 

Finally, as an aid to simplifying the design procedure, it is 
sometimes assumed that vertical and pitch dynamics also decouple. This 
is a much cruder approximation and is retained only to make certain 
design decisions and not for the design evaluation. 

3.2.3 Optimization Problem 

The optimal controller for the tracking problem is presented in 
this section. Since the system is linear, with additive gaussian white 
plant noise, and the performance index can be expressed as a quadratic 
cost function defined over an infinite time interval; the steady state 
linear quadratic gaussian (SSLQG) regulator is the optimal controller. 

The SSLQG optimal controler problem and solution is given in ref- 
erences (23) and (24). It can be stated as follows: 

Find a feedback control law jj(t) = ju(x(t)) to minimize 
the quadratic cost functional, 

T 

J(u) = E { lim | J* (£ T (t)Q£(t) + ii T (t)Ru(t)dt } 

T-*» 1 u 

(3.2. 3.1) 


(3. 2. 3. 2) 


subject to the constraints , 
x(t) = Ax(t) + B_u(t) + B Jf £(t) 
£(t) = Cx(t) 


(3. 2. 3. 3) 
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where £A, B). is atabilizable , 

(A, C) is detectable , 

P is a gaussian white noise process, 

and E { • } is the expectation operator. 

T 

Q = Q > 0 is the per unit time quadratic cost associated with 

the response of y(t) , 

T 

R =* R > 0 is the per unit time quadratic cost of the control 
u(t). 

The form of the feedback control ii(x(t)) is also given in the above 
references. u.(x(t)) follows: 

The cost functional J is minimized by, 

u(x(t)) * Gx(t) - -R -1 B T P x(t) (3. 2. 3. 4) 

where P is the unique symmetric positive semidefinite solution of 
the matrix algebraic Riccati equation (MARE) , 

0 - PA + A T P - PBR** 1 B T P + Q (3. 2. 3.5) 

Note that for given system state and response equations 
(3. 2. 3. 2) and (3. 2. 3. 3), an SSLQG regulator design is uniquely deter- 
mined by the response and control vector quadratic cost weighting 
matrices Q and R. Also note that the control problem formulated above 
is a special case of the more general one where A, B, B^, C, Q, and R 
are time varying. 

The implementation of the feedback gains results in a closed loop 
system whose state equation is written: 
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x = (A - BR _1 B T P)x + (3. 2. 3.6) 

or using the notation: 

A a =* A - BG - A - BR -1 B T P (3. 2. 3. 7) 

the state equation becomes: 

i - A a x + B^£ 

In the overall system of the aircraft, ship and waves; only the 
aircraft is controllable. Therefore, the closed loop system matrix, A 
has the form shown in Table (3. 2. 3.1), where: 

A^ is the aircraft open loop system matrix, 

Ag is the ship system matrix, 

Aq^ is the ocean wave system matrix, 
is the aircraft input matrix, 

G^ is the matrix by which the aircraft states are fed back 
to the aircraft controls, 

Gg is the ship state vector feedforward matrix to the aircraft, 
Gq W is the ocean wave state vector feedforward matrix to the 
aircraft* 

Figure (3. 2. 3.1) is a block diagram of the hover controller. In the 
figure, the ocean wave and ship models are combined for compactness. A^., 
B^g and G^g represent the combined ocean wave and ship state vector feed- 
forward matrix. 

For the hover control problem having only vertical and pitch 
dynamics, the following correspondences are also made to the general 
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OVERALL CLOSED LOOP SYSTEM MATRIX 
A a - [A - BG] 

0 ‘ C °A I G S I C 0» 3 



a a + b a g a 

1 

, b a g s 

! 

I b a g ow 

a cl ■ 

0 

1 

A <? 

■1 

I ^s 


0 

1 

0 

i 

1 A ow 


TABLE (3. 2. 3.1) 
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control problem stated above. _p Is the zero mean, gaussian, white noise 

driving the wave model, u is the control vector driving the aircraft. 

Note, a has two elements, altitude control and pitch control. In the 

quadratic weighting matrix R, these are weighted by r^ z and r^g . The 

A A 

response for ship tracking is chosen as shown: 


z(aircraft) - z(ship deck) 

=3 

Z A-S 

0 (aircraft) - 9 (ship deck) 


9 a-s 


( 3 . 2 . 3 . 8 ) 


Therefore, the controlled variables are the vertical position and pitch 

orientation tracking errors of the aircraft with respect to the ship. 

The corresponding diagonal elements of Q are q_ and q fl . The 

Z A-S y A-S 

weighting q_ causes the vertical position of the aircraft to match 
A-S 

that of the ship. A constant non-zero separation, however, does not 
change the system dynamics, so g may be viewed as the deviation of 
the aircraft from its desired position with respect to the ship. q fl 

A-S 

likewise causes the aircraft to match the ship's pitch attitude. 


3.3 Controller Design and' Evaluation Method 
3.3.1 Design Method 

This section discusses how the assumptions and theory are employed 
in the design procedure. The design procedure begins with the design 
of a basic hover controller for the case of a stationary ship deck. That 
is, for the case when the ship has no motion. This design is developed in 
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two steps. The first step is to gat a controller for the pitch dynamics 
only. The second step is to use this design for pitch, and determine the 
additional controller structure needed to get a satisfactory controller 
for coupled pitch-heave dynamics. The basic design's quadratic weightings 
are then used to handle the ship motion tracking problem. The aircraft 
model is augmented by the ship and wave models, and the. aircraft altitude 
and pitch quadratic cost weightings are transfered to the aircraft-ship 
separation and orientation outputs. By varying the two output weightings 
in unison, controller designs are obtained that reflect the relative 
importance of tracking performance compared to control effort cost. 

Design method details follow. In the basic design, both pitch and 
vertical control are weighted by unity. For complete design freedom, 
these could be varied independently. However, since the pitch and verti- 
cal dynamics of the aircraft are approximately uncoupled, the pitch-heave 
control problem can be separated into a vertical control problem and a 
pitch control problem. In each of these, the important quantity is the 
relative weighting of output to input. Hence, the choice of unity con- 
trol weightings. Note that in particular, the vertical dynamics do not 
affect the pitch dynamics, so the results of step two do not alter those 
of step one. 

The two step design method generates a basic controller that is 
completely determined by two output and two control weightings. The 
basic design is for hover control, that maintains the aircraft in a 
horizontal attitude over a fixed ship deck at a constant vertical sepa- 
ration. The commanded separation is defined to be zero for convenience. 
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A suitable basic pitch control design is one having a response band- 
width typical of piloted fighter aircraft. Reference ( 9 ) gives pilot 
ratings versus aircraft handling qualities. Bandwidths suitable for 
manual control were chosen, because aircraft behaving significantly 
different under automatic and manual control, or unlike conventional 
aircraft, may receive poor pilot acceptance. In fact, similar manual/ 
automatic behavior may be a safety consideration, required to prevent 
undesirable transients in situations where the pilot is forced to switch 
from automatic to manual control. Of course, for good ship motion track- 
ing, the bandwidths must also be at least as large as the dominant ship 
motion frequencies. 

For making bandwidth ’estimates , the Butterworth pole pattern is 
useful. An explanation follows. Pitch and heave dynamics are two 

approximately uncoupled, third order, single input-single output sytems 
(pitch control input and pitch output and altitude control input and 
altitude output). (As usual, all states are assumed to be available for 
feedback control.) Neither has transfer function zeroes, and so the 
SSLQ regulator produces third order Butterworth pole patterns for each. 
These patterns are generated by varying the relative quadratic cost 
weighting magnitudes of the input and output. As the output weighting 
is increased, the poles approach a common magnitude and move radially 
outward from the origin along asymptotes. A detailed explanation can be 
found in reference ( 23 ) • Since the poles of such Butterworth patterns 
all have the same magnitude, this magnitude is used as a convenient 
measure of system bandwidth. In addition, the Butterworth pattern is 
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useful for design because it indicates that by increasing the output 
weighting, while holding the input weighting constant, the closed loop 
system disturbance rejection bandwidth can be made as wide as desired. 

To adapt the basic design to track a moving ship deck, the system 
is first augmented to include the ship and wave states of interest. For 
the ship, pitch and heave states are required, and for the waves a com- 
plete model with the pitch and heave outputs is needed. Secondly, the 
penalty weightings are transfered to the aircraft to ship relative pitch 

and vertical separation. They become q 0 and q„ . The weightings 

°A-S a-S 

determined for the basic design now only indicate estimates of the ones 
required. However, the determination of these basic weightings provides 
useful experience in estimating the magnitude of weighting changes in 
the quadratic cost functional and the corresponding performance changes. 
The final choice of q Q and q_ is made by scaling them both simul- 

9 a-s z a-s 

taneously. Since these are the only non-zero- elements of the Q matrix, 
this is just a scaling of Q. Scaling Q to get a final design is an 
arbitrary, but convenient, choice, since it reduces the design parameters 
down to one. In addition, by varying all of Q simultaneously, the rela- 
tive importance of pitch and heave performance is maintained. 

Having described the design method, it is appropriate to discuss 
how designs are evaluated. This is the topic of the next section. 

3.3.2 Performance Evaluation Method 

Quantities like tracking error, control effort, and bandwidth are 
of interest in deciding the acceptability of a tracking controller. 
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Several of these can be estimated from the steady state covariance 
equation, while bandwidth can be obtained from the system pole locations. 

The list of quantities calculated for controller performance evalu- 
ation follows: 

(1) rms aircraft to ship relative vertical separation: 

a„ feet 
Z A-S 

(2) rms aircraft to ship relative vertical velocity: 

feet/second 

A-S 


(3) rms aircraft to ship relative pitch orientation: 


degrees 


'A-S 


(4) rms aircraft to ship relative pitch rate: 

0 ’s degrees/second 
0 A-S 

(5) rms aircraft vertical control actuator authority: 

percent of nominal 


(6) rms aircraft vertical control rate of change: 

percent of nominal/second 

(7) vertical loop bandwidth: radians /second 

The first six of these quantities are calculated from the steady 
state covariance equation for the outputs and states. In evaluation, 
three root-mean square (3rms) values are used, since they are rarely 
exceeded. The last quantity, vertical loop bandwidth is the magnitude 
of the closed loop system eigenvalues associated with the vertical 
dynamics. Recall that this bandwidth estimate is by virtue of the 
Butterworth pole configuration discussed in the previous section. The 
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pitch loop bandwidth could be determined similarly. However, it is found 
that .the pitch loop requirements of this system are relatively easy to 
satisfy compared to the vertical loop ones. As a result, only orientation 
and rate behavior are recorded for the pitch axis. 

The steady state roo t-mean-square deviations of the outputs and 
states are determined by the following method ( 23 , 25 ) . Let be the 
vector of quantities whose rms values are desired. Write in terms of 
the state vector x as: 

Z - Cx (3. 3. 2.1) 

The the covarinace matrix of £ is given by: 

cov(£) = CPC T (3. 3. 2. 2) 

where 

cov(x) * P (3. 3. 2. 3) 

In turn, P. is the steady state solution of the covariance equation: 

0 = P = A CI P + Pa£ l + BVB T (3. 3. 2. 4) 

where A^ is the closed loop system matrix, 3 the noise input matrix, and 
V the covariance intensity of the driving, zero mean, gaussian, white 
noise £ of the closed loop system: 

x = A CL— + (3. 3. 2. 5) 

The lrms, 1-sigma, or 1 standard deviation values of are just the square 
roots of the diagonal elements of cov(y) . 1-sigma values, however, are 
not themselves indicative of the maximum excursions of the elements. 

From probability theory, it is known, that for a gaussian random variable. 
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the percentage of outcomes within 1, 2 and 3 standard deviations of the 
mean is 68.3%, 95.4% and 99.7% respectively . (computed using reference 
( 26)). For ergodic stochastic processes, "percentage of the outcomes" 
can be substituted with "percentage of the time". For the plots presented 
in the design evaluation, 3rms values are chosen. These values are ex- 
ceeded only 0.3% of the time, and are therefore conservative estimates 
of maximum or ultimate levels for design specifications. 

3.4 Controller Design and Evaluation 
3.4.1 Introduction 

The following sections describe the hover controller design and 
evaluation. The first section discusses the basic controller for hover 
at a fixed point. As described earlier, this design is developed in 
two steps. The first step is the choice of a pitch dynamics controller. 

In the second step, vertical dynamics are included and a correspondingly 
augmented controller is chosen. The basic controller is converted and 
used as an initial ship motion tracking controller design in the subse- 
quent section. That section reports on the tradeoff between tracking 
performance and control cost obtained from scaling the output weighting 
matrix of the quadratic performance index. Finally, a section is in- 
cluded on the effect of ship modelling accuracy. There, the performance 
results for a low order ship motion model are compared to the results 
obtained for the standard wave-ship model. 
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3.4.2 Quadratic Weights for the Basic Design 

The two step scheme outlined in the design method section was used 
to choose a basic design for the hover control regulator. 

I 

Step one is the design of a pitch controller. This is designed as a 

single input controller with full state feedback by weighting the pitch 

control input and the pitch attitude output in the quadratic performance 

index. For each ratio of these weights, a set of feedback gains and 

closed loop pole locations is determined. The Butterworth pole pattern 

for this system is shown in Figure (3. 4. 2.1). The input weighting is 

unity as explained in the design method section, and the weight q s is 

A 

10.0. This gives closed loop pole location magnitudes and bandwidths of 

f 

about 3 radians/second, which is typical of acceptable piloted aircraft 
responses ( 9 ) • 

Step two is the inclusion of vertical control. The system is now 

sixth order, with two inputs and two weighted outputs. However, the 

coupling between pitch and heave is not strong, so the choice of q_ , 

Z A 

the vertical position state weighting, does not influence the closed 
loop eigenvalues of the pitch dynamics. The pole location pattern for 
the vertical dynamics is shown in Figure (3. 4. 2. 2). The vertical posi- 
tion weighting, q„ , is chosen to be 0.01, giving a closed loop vertical 
Z A 

response bandwidth of about 1.5 radians /second (Reference ( 7 ) suggests 
similar attitude and flight path dynamic bandwidths at 2.0 and 1.25 
radians /second respectively.) The fact that vertical bandwidth is less 
than the pitch bandwidth reflects physical constraints. Attitude control 
power of the VTOL aircraft is obtained by rerouting power between the 
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fans, without requiring any change in the overall engine output power. 

This can be done relatively quickly. However, vertical control requires 
varying the engine output power, which is a slower process. In addition, 
vertical thrust is limited by the aircraft weight and maximum available 
engine thrust. This last constraint is expressed as a thrust to weight 
ratio, T/W. For stationary, undisturbed hover T/W is 1.0. Typically, 
experimental VTOL vehicles have a maximum T/W capability on the order of 
1.05 to 1.3 (27 , 28 ) . 

The two weightings chosen- above specify the initial point or basic 
design for the tracking controller. The latter has aircraft-ship rela- 
tive pitch and vertical separation outputs as controlled responses. Thus, 
the initial quadratic weightings for the tracking hover controller are 

q Q =10.0 and q„ = 0.01. 

0 A-S Z A-S 

3.4.3 Tracking Controller Performance 

This section surveys the performance characteristics of the ship 
motion tracking hover controller. The basic design of the last section 
is the starting point for a study of the relation between tracking per- 
formance and control cost as quadratic ratios are varied. The covariance 
analysis and bandwidth estimation scheme described earlier shows how- 
tracking errors go to zero and control effort becomes large-,- while the 
bandwidth and the weightings ratio increase. 

These results are illustrated in the plots of Figures (3. 4. 3.1), 

(3. 4. 3. 2) and (3.4.3. 3). Each of these plots has aircraft vertical 
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control actuator authority magnitude as the abscissa, and one of the other 
six quantities listed in section (3,3,2) as the ordinate. Control actu- 
ator authority is presented as aircraft thrust to weight ratio. This is 
defined as 1.0 (the thrust to weight for stationary, undisturbed hover), 

plus the 3rms actuator level (in fractions of the amount for stationary 

• • 

hover). The ordinate quantities , z^_ s , ®k.-S ’ ant * vert * ca l 

thrust rate of change are plotted as 3rms values. The vertical control 
bandwidth is taken as the vertical actuator eigenvalue magnitude in 
radians /second. The index along each curve is q, the ratio q equal to 1.0 
gives the basic design. The output quadratic weighting matrix is given: 

Q * DIAG [ 0, q fl , 0, q^ , 0, ...» 0 ] (3. 4. 3.1) 

9 a-s z a-s 

where 

q Q * q • 10.0 (3. 4. 3. 2) 

9 A-S 

q = q • 0.01 (3. 4. 3. 3) 

Z A-S 

The control input quadratic weighting matrix is given: 

R = DIAG [ r. a , r. ] = DIAG[ 1.0, 1.0 ] (3. 4. 3. 4) 

60 A-S 6z A-S 

The first four of the plots indicate that as thrust to weight (T/W) 
is increased, the aircraft hover tracking performance improves. This 
trend corresponds to increasing the output weighting matrix magnitude 
q: the more heavily the position and orientation errors are penalized 

with respect to control, the smaller these errors are kept, and the 


more control and thus actuator effort is used. 
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Figure (3. 4. 3. 3) indicates that at the same time tracking errors 
decrease and thrust to weight increases, the rate of change of thrust to 
weight (or actuator rates) increases. Similarly, the bandwidths of the 
aircraft dynamics increase. 

The T/W ■ 1.0 point on the curves represents the condition when the 
aircraft is hovering at a fixed altitude. Separation errors are not 
weighted at all at this point, so vertical separation, relative velocity, 
pitch, and pitch rate are due only to ship motion. Thus, for example, 
the 3rms ship motion is 12.5 feet (above and below the mean). 

The vertical control bandwidth curve has an interesting shape. It 
is quite flat near T/W » 1.1. In this region, the bandwidth of the 
aircraft vertical motion is between 0.5 and 1.0 radians/second. To ex- 
plain this, note that the magnitude of the closed loop pole locations 
increase with the output weighting. On the other hand, T/W only increases 
at frequencies where ship motion has significant power. Therefore, band- 
width increases at low frequencies, but T/W does not. At ship motion 
frequencies, T/W increases significantly causing the curve to flatten. 

Then at high frequencies, where ship motions are small, T/W increases 
very slowly; however, bandwidth continues to grow and the curve rises 
sharply again. 

To demonstrate tracking improvement due to increased thrust to weight, 
simulations are shown in figures (3.4. 3. 4) and (3. 4. 3. 5) for two different 
output weighting magnitudes. Table (3. 4. 3.1) gives the parameters of 
these simulations. In the figures, the aircraft and ship motions are 
superimposed, so the difference represents tracking error as usual. (In 



Design Parameters for the VTOL Tracking Simulations in Figures 3. 4. 3. A and 3. 4. 3. 5 


Simulation 

Output Weighting 

Control Weighting 

Thrust to Weight 

Tracking error 


% 

A-S 

m 

r 60 

r «sz 

T/W 

(1+3 rms) 

Z A-S 

(3 rms, feet) 

1 

0.1 

0.0001 

1.0 

1.0 

1.08 

6.9 

2 

1.0 

0.001 

1.0 

1.0 

1.16 

2.4 


Table 3. 4. 3.1 
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Tracking Simulation 1: Thrust to Weight = 1.08 

Figure 3. 4. 3.4 
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a physical realization, a non-zero separation would be included, which 
would be slowly decreased to zero to achieve a landing.) As the table 
shows, the factor of 10 increase in the output weighting between these 
two simulations leads to a thrust to weight increase of a factor of 2, 
from 8% to 16% control authority. At the same time, the corresponding 
vertical tracking error decreases by a factor of about 3. The simulation 
plots clearly illustrate this tracking improvement. 

3.4.4 Sensitivity to Ship Modelling Accuracy 

This section presents evidence that the accuracy of ship motion 
modelling significantly affects the performance capabilities of the 
tracking controller. 

The performance curves obtained to this point have been for a ship- 
wave model whose longitudinal dynamics are tenth order, and approximate 
ship responses to waves having a Neumann ocean rwave' spectrum. However, 
before this model was developed, a less sophisticated fourth order model 
was used. Both the tenth order "standard model" and the fourth order 
"low order model" are described in Chapter Two. 

Figure (3. 4. 4.1) gives the vertical separation error versus thrust 
to weight for tracking the low order model. The first difference to 
be noted is that at T/W = 1.0, the separation, again due to ship motion 
alone, is almost 17 feet (3rms) rather than the 12.5 feet (3rms) for the 
standard model. This is because pitch and heave motions in the low order 
model are independent, and so their effect on the landing deck motion 
variance is additive. In the standard ship-wave model, these motions 
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are correlated in a way that yields less motion at the landing deck. The 
other difference is the rate at which tracking error is decreased with 
increased thrust to weight. This rate is very slow for the low order 
model. The reason for this behavior is that the spectrum of ship motions 
for the simple model is not nearly so peaked as for the sophisticated 
model, and so there is a much wider spread of frequencies over which 
ship motions are significant. Thus, the low order model requires higher 
values of T/W to arrive at the same level of tracking performance as the 
standard model. A sketch of the magnitude spectra of the two models was 
given in Figure (2. 6. 6.1). 

Figure (3. 4. 4. 2) shows the vertical control bandwidth for the low 
order model. Compared to the corresponding figure (3. 4. 3. 3) for the 
standard model, the vertical control bandwidth increases much more slowly 
with thrust to weight, indicating the wide range of frequencies over 
which the low order model has significant motion energy. In addition, 
the peak motions now occur near one radian/second. They also have larger 
amplitude requiring a larger increase in aircraft T/W to track them. 
Although plots are not given, pitch tracking performance shows similar 
trends . 

These comparisons of the ship tracking performance capabilities 
of the hover controller demonstrate the importance of accurate ship 
motion modelling for meaningful assessments of aircraft hovering per- 
formance, and associated control system bandwidth and thrust to weight 
requirements . 
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3.4.5 Conclusions 

The goal of the hover controller was to make the aircraft track the 
ship deck motion as closely as possible, using only a small amount of 
vertical actuator control authority. The results obtained are perfor- 
mance curves illustrating the trade-off between tracking errors and 
required aircraft thrust to weight ratio, and the relationship between 
closed loop system bandwidth for vertical dynamics and the thrust to 
weight ratio. These curves are indicative of the "best" performance 
achievable by the VTOL aircraft. 

Although only vertical and pitch dynamics are considered, the design 
method and analysis technique could be used for generating a complete 
full dynamics tracking controller. 

Modelling accuracy was found to be very important in obtaining 
meaningful results. A low order ship model having a significantly 
wider power spectrum than the more accurate standard model, required a 
controller with much higher actuator authority than the standard one. 

The results of this chapter motivate the landing controller design 
in Chapter Four. To meet a goal of reducing touchdown velocities to 
about one foot per second or less, the hover controller requires a 
thrust to weight ratio greater than 1.1 (figure (3. 4. 3.1)). However, 
the maximum ratio of current VTOL aircraft is only about 1.3 (27 , 28 ) . 
Thus, the control effort used by the hover controller is a large frac- 
tion of the available authority. This motivates the search for an al- 
ternative design to achieve good landing performance at lower actuator 
control levels. Chapter Four discusses one such alternative design. 
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CHAPTER 4 
LANDING CONTROL 


In this chapter, a landing controller design for the VTOL aircraft 
is developed. The purpose of the controller is to land the aircraft on 
the ship deck beginning with the aircraft hovering at some initial alti- 
tude. The design of this chapter is motivated by the performance of 
the hover controller in the previous chapter. There it was found that 
simple ship deck tracking requires potentially high levels of control 
authority to achieve low relative velocities suitable for landings. The 
objective of the landing controller in this chapter is to reduce this 
thrust to weight requirement. 

This chapter is divided into four main parts. The introduction 
presents the basic concept of the proposed landing controller. Part two 
states an optimal control problem whose theoretical solution implements 
the concept. Only a suboptimal solution has been obtained for this 
problem. These are developed in part three, the analysis and design. 
Finally, the last part gives an evaluation of the resulting controller. 

The goal of the new controller is to reduce the control actuator 
levels required to land, below those required by the tracking controller. 
To achieve ‘this. 'goal, the new controller exploits the predictability 
of ship motion and lands the aircraft at a peak of vertical ship deck 
motion, as sketched in Figure (4.1.1). A simple calculation shows how 
l.uw the required actuator level could be in a situation where the peak is 
predicted perfectly. The calculation is just the control effort required 
by a bang-bang controller to move a mass between two paints in a given time 
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with zero initial and final velocities. Gravity is not a factor in the 

calculation, since the aircraft model is linearized about stationary 

hover. The calculation goes as follows. Denote the initial and final 

times of the trajectory by t and t , respectively. During the first 

X r 

half of the trajectory, from time t_ to i(t - t ), the acceleration of 

the mass is +a. At the end of the interval, its position will be 
11 2 

-^a(^(t F - t^)) . During the last half of the trajectory, the acceleration 
will be -a, and the mass will have moved an additional distance equal to 
the first. When the mass comes to rest, the total change in position 
AS will be: 

AS - \ a (t p - t x ) 2 (4.1.1) 

The acceleration level to move the mass the distance AS is therefore: 

a = 4 AS/ (tj, - tj) 2 (4.1.2) 

For tp - tj » 10 seconds (the natural period of the ship is roughly 9.6 
seconds), and AS = 8 feet (for hover at 13 feet above mean sea height 
and descending to a peak 5 feet above the mean height) , the acceleration 
magnitude is |a| =0.32 feet/second . This corresponds to a control 

actuator level of about one percent of the nominal stationary hover 
requirement. The hover controller of the last chapter required an actu- 
ator authority roughly an order of magnitude larger than this value. 

In addition to the potentially reduced thrust requirements, other 
benefits may be realized by this landing control strategy over the hover 
controller. The first advantage is that ship and wave states are not 
required aboard the aircraft. The only information needed about the 
ship is the prediction of the time and position at which the landing pad 
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will reach a crest. Compared to the ship state vector, this is less 
complex information to transmit. Furthermore, it may be possible to 
transmit the information at a lower rate than would be necessary for 
tracking control. The required update rate remains to be determined 
however, before the magnitude of this benefit can be ascertained. 


4.2 Control Problem 


4.2.1 Optimization Problem Statement 


The cause of high thrust to weight ratios in the aircraft ship 
motion tracking controller is the presence of significant power at high 
frequencies in the motion spectrum of the ship. The goal of the new 
controller is to eliminate this effect. Ideally this controller com- 
mands the aircraft from stationary hover at some initial height along 
a trajectory of low vertical actuator control levels, towards a height 
and time when the ship deck motion will reach' a crest, ending with zero 
relative aircraft-ship deck separation and velocity. This ideal behavior 
can be expressed in terms of the following quadratic optimal control 
problem. 

Find a feedback control law _u(t) =* £(x(t)) to minimize 
the quadratic cost functional 

J(u) = E{ £ T (t F ) P F £(t F ) + / F /(tOQCtO^Ct) + u T (t)R(t)u(t) dt } 

1 (4. 2. 1.1) 


subject to the constraints 


x(t) = Ax(t) + Bu(t) + B N £(t); xCtj) = Xq 


(4. 2. 1.2) 
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X(t) - C^xCt) 

where the terminal time is defined by 

t F » minimum {argument [minimum C 2 x(t)]} 
T t > t- 


(4.2 .1.3) 


(4. 2. 1.4) 


P_ =* P > 0 is the quadratic cost associated with the 

r r 

end time response £(tj.) , Q(t) = Q^(t) > 0 and R(t) = R^(t) > 0 
are the per unit time costs associated with the response _jr(t) and 
control jj(t) respectively; and (A,B) is stabilizable and (A,C) is 
detectable. 

In the quadratic cost functional P^, Q(t)and R(t) should be chosen so 
the resulting aircraft trajectories use low vertical control effort 
and have small terminal errors £(tp) . 

This stochastic control problem is the standard finite terminal 
time problem, except that the performance index J is a random variable 
depending upon the ship motion, and thus on £(t) for tg <_ t <_ t^ through 
the relation of equation (4. 2. 1.4). While a formal solution to this 
optimal control problem with random end time is not known, the following 
heuristic design provides a reasonable approximate solution. 


4.2.2 An Approximate Solution 

To take advantage of existing results for linear quadratic control 
problems, the solution of the optimization problem in Section C.4.2,1) was 
approximated by treating the best currently available estimates of 
terminal time and corresponding deck position as fixed deterministic 
parameters; An explanation of the resulting algorithm follows. 
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The aircraft is initially assumed to be at stationary hover a safe 
distance above the ship landing deck. Ship motion is predicted ahead 
to the next crest. The time and height of this crest are compared with 
the current time and aircraft altitude to determine if a landing on the 
crest is possible. If not, the aircraft is commanded to remain at its 
original position. Otherwise, a landing is initiated. Using the pre- 
dicted crest time and height, a nominal aircraft trajectory to the end- 
point is computed via the standard time varying linear quadratic regul- 
ator problem. This trajectory is then followed by the aircraft using 
a steady state linear quadratic gaussian (SSLQG) regulator to control 
trajectory deviations. At certain intervals in time, new predictions 
of the ship motion peak are made and the nominal trajectory updated. 

All the gains used in this design can be calculated prior to the 
landing. The trajectory following SSLQG gains used for tracking the 
desired trajectory are, of course, constant and precomputable. Gains 
for the nominal trajectory itself are precomputable and independent 
of the vertical height through which the aircraft must move, the initial 
vertical velocity of the aircraft, and the time required for the man- 
euver. They are only functions of time-to-go and thus, can be stored 
as a single time varying set of gains. During a landing approach, these 
stored nominal gains are used to compute a new nominal trajectory each 
time an updated estimate of the end point and time are available. 
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4.3 Analysis and Design 

4.3.1 Assumptions 

Here the additional assumptions beyond those posed for the modelling 
and the hover controller are stated. All the previous assumptions leading 
to the linearized model of the aircraft-ship system, and the design of 
the longitudinal hover controller in the previous chapters are required. 
For emphasis, it should be remembered that the ability to decouple air- 
craft pitch and heave dynamics was used and aircraft states were con- 
sidered to be perfectly known. In addition, predictions of ship landing 
pad motion extrema are now required. In the simulations used for eval- 
uation of the controller, peaks are identified via a state space model, 
but in fact the method of prediction is not important. Ship and wave 
state measurements are not required as key elements of the algorithm. 

4.3.2 Touchdown Prediction Analysis 

An approximate analysis of touchdown errors is made in this section. 
The results are given in Table (4. 3. 2.1). The expected velocity errors 
are used in section (4.3.6) to obtain a bound on the nominal aircraft 
landing trajectory update frequency, and the expected time errors are 
compared with the simulation results of Section (4.4.1). 

A touchdown occurs when the relative aircraft-ship velocity is 
greater than zero at the first moment when their separation is zero: 
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EXPECTED TOUCHDOWN ERRORS 


PREDICTION TIME 
C F 

° A W 

v s (t F > 

EXPECTED 

TOUCHDOWN 

VELOCITY 

EXPECTED 
TOUCHDOWN 
TIME ERROR 

seconds 

feet 

feet/second 

feet/ second 

seconds 

1.0 

0.05 

0.14 

0.45 

0.18 

2.0 

0.17 

0.20 

0.80 

0.32 


Table 4.3. 2.1 



> 0 
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(4.3. 2.1) 

(4.3.2. 2) 


Z A-S ( V 
Z A-S^ V 


0 


A touchdown is characterized by its time, the vertical position of 
the ship deck, and the relative aircraft-ship vertical velocity. The 
landing controller is designed to generate a near zero velocity touch- 
down between the aircraft and ship by commanding the aircraft along a 
trajectory, which terminates with zero velocity and meets the ship deck 
at its predicted crest height and time. However, the actual touchdown 
has different values due to uncertainties in the aircraft and ship 
trajectories. In the present model, the aircraft motion is deterministic 
once it has been given a commanded trajectory, so all the uncertainties 
are presently due to ship motion only. 

The relative separation and velocity of the aircraft and ship can 
be written as the sum of their deterministic and stochastic components 
as follows: 


(t) 

= E {z A _ s (t)> 

+ ^ k _ s M 

(4.3. 

.2.3) 

(t) 

= E{ z A _ s (t)} 

+ A Vs (t) 

(4.3. 

.2.4) 


Equation (4. 3. 2. 3) is illustrated in Figure (4. 3. 2.1). 

Rewriting these equations with the bar notation for mean, noting 
that the uncertainties are due only to the ship (Az^ g (t) = - Az g (t), 
Az A _ s (t) = -Az s (t)), and then dropping the subscripts yields: 


-z A _ s (t) = z(t) + A z(t) 


(4. 3. 2. 5) 


■ Z A-S (t) 


z(t) + Az(t) 


(4. 3. 2. 6) 



VERTICAL POSITION 
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Deterministic and Stochastic Components of a Landing Trajectory 
Figure (4. 3. 2.1) 
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«■> • 

z(t) and z(t) are approximated by 

2(t) = 0(t - t F ) 2 

z(t) ■ . 2a(t - ty) 


(A. 3. 2. 7) 


( 4 . 3 . 2 . 8 ) 


where a » 1.259. This approximates to second order a cosine wave that 
peaks at the predicted endtime t^., and has the same typical period (9.6 
seconds) and rms amplitude (4.2 feet) as ship deck motion. For the 
ideal case where the aircraft has reached the end of its trajectory and 
is nearly stationary, the approximate z(t) and z(t) represent the total 
expected aircraft-ship separation and relative velocity. As is clear 
from Figure (4. 3. 2.1), actual touchdowns will occur early and at slightly 
negative positions. The deviations existing at touchdown time t^ can be 

approximated by the values of t that would occur "if the aircraft and 

IT 

ship could pass through each other after touchdown". This is illustrated 

by the dotted line. If for the interval between t and t , the velocity 

T r 

deviation Az(t) is assumed to be constant, then a reasonable approximation 
for the touchdown deviations are given in terms of deviations at t by: 

r 

Az(t F ) = Az(t T ) + A£(t F )(t F - t ] .) (4. 3. 2. 9) 

Az(t F ) = Az(t T ) (4.3.2.10) 


The touchdown time can be found by substituting equation (4.3.2.10) 
into (4. 3. 2. 9), then substituting the result and equation (4. 3. 2. 7) into 
(4. 3. 2. 5) with z^ ^(t^) equal to zero: 


C F " C T 3 (Az Ct-p) + (Az 2 (t F ) - 4aAz(t F )) 1/2 /2a (4.3.2.11) 
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The touchdown velocity is obtained by substituting this last expres- 
sion along with equations (4. 2. 3. 8) and (4.3.2.10) into (4. 3. 2. 6): 

Vs ( V = ~Az(t F ) - 2o(t T - t ? ) (4.3.2.12) 

- (Az(t F ) - 4ctAz(t F )) 1/2 (4.3.2.13) 

To second order in Az(t F ) this is: 

*z A _ s (t T ) = (~4ctAz(t F ) ) 1/2 + |(-oAz(t F ))“ 1/2 Az 2 (t F ) (4.3.2.14) 


The expected touchdown velocity is found by evaluating: 


E { VsM Vs'V =• °> Vs ( v - 0 > (4.3.2.15) 

The conditions on z^_^ ( t^) and z^_^ (t^) in this expectation are satis- 
fied by all trajectories such that Az(t_) < 0 (plus some other trajec- 

r 

toreis, which are neglected). Hence, 


E {^(tp | Vs ( V =• °> Vs^ = 0 } 


= E {(-4o4z(e p )) 1/J + |(-cc4z(t F ))' 1/2 4z(t F ) |Az(t F ) < 0 } 

(4.3.2.16) 

The expectation is over Az(t F ) < 0 and all Az(t F ), which are gaussian 
random variables with zero means and variances depending on t F » where 
t F is measured with respect to the time when the ship state vector was 
last updated. Statistics for Az(t_) and Az(t_) are computed at the end 

r r 

of this section. They are approximately independent with standard 

deviations a. and a. • . . known as functions of t„. 

Az s (t F ) Az s (t F ) F 
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The resulting equation for the touchdown velocity can be written 
as : 


E { E 

Az ( ty) Az(t F ) 


{ (-4aAz(t F )) 1/2 + |(-aAz(t i ,))" 1/2 Az 2 (t F )|Az(t F ) < 0 }} 

(4.3.2.17) 


• . • 

but the expectation of Az (tj.) with respect to Az(t ? ) is just a ^z(t )’ 

F 

and the probability density for Az(tj,) over non-positive values is 

*Az|Az < 0 (Z I Z < °> - exp{ ~2 (Z/0 Az(t I ) >21 

F F (4.3.2.18) 

Therefore the expectation reduces to: 


E { (-4aAz(t 1? )) 1/2 + ^(-aA2(t F ))" 1/2 a 2 * (t ^ ) | AzO^) < 0 } 

- ^ d/0 42(tj) > V° [ (-4ctZ) 1/2 + i(^z)' 1/2 a 4 - (tr) ] 

exp{ - |(Z/a Az(t } ) 2 } dZ (4.3.2.19) 

F 


The two terms of the integral are given in reference ( 26 ) . Evaluating 
the constants gives: 


E{ VsMWtP 5 °- ‘w'V • 0 } 


» 1M5C 1 / 2 ,, . + 0.3832(0. .. ,)' 1/2 0!-,. . 

Az(t F ) Az(t p r Az(t F ) 


(4.3.2.20) 


Numerical values for this expression are given in Table (4.3. 2.1). 
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The expectation of touchdown time error t^, - t^, can be found by 
noting that the expectation of Az(tj.) in equation (4.3.2.11) is zero. 
The remaining term in (4.3.2.11) is the expectation of the touchdown 
velocity divided by 2a. Hence, 


E { t_ - t_ z. _ (t_) > 0, z. _(t_) =» 0 } 

F T A-S T A-S T 


= (X^2a)E{ z A _S (t T ) l Z A-S (t T ) > °’ 2 A-S (t T ) = 0 } (4.3.2.21) 

Table (4.2. 3.1) gives these expectations as well. 

An expression for the expected ship deck position is not required 

since in this analysis all collisions occur at z g (t^,) equal to zero. 

In actual touchdowns, however, the expected ship deck position is often 

above the predicted height at touchdown since the aircraft trajectory 

is designed such that it does not dip below the predicted end point. 

The supporting calculation of the variances of Az_(t_) and Az_(t x .) 

or or 

follows . 


The growth of ship motion prediction uncertainties (a^ z [\ z ( t )^ 

S S 

is calculated by solving the time varying covariance equation for the 
ship deck vertical position and velocity. The state space ship model 
was used to determine the quantities required in order to solve the co- 
variance equation 1 '.'' The solution P(.t) is the state covariance matrix for 
the ship and wave states at time t. The initial condition for the solu- 
tion P(0) is equal to the zero matrix. This indicates that all the 
states are perfectly known at time t= 0. The matrix differential equa- 
tion (given earlier in equation (3. 3. 2. 4)) is: 
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P(t) = AP(t) + P(t)A T + BVB T (4.3.2.22) 

with initial condition: 

P(0) » 0 (4.3.2.23) 

A, B and V are respectively system, noise input, and noise intensity 
matrices for the wave-ship model. The output covariance matrix was 
then found from: 

Cov(y (t) ) = C(t)P(t)C T (t) (4.2.3.24) 

For this analysis, the outputs were ship landing pad vertical position 
and velocity. The results are shown in Figures (4. 3. 2. 2), (4.3.2. 3) 
and (4. 3. 2. 4). These give rms position and velocity, and the correlation 
of position and velocity. 

Since the output y(t) is the result of a linear system driven by 
a gaussian process, it too is gaussian, and so 68% of all outcomes are 
within 1 rms (one standard deviation) of the mean, 95% within 2 rms, 
and over 99% within 3 rms. (Recall this result was also used in Section 
3.3.2).) Thus, in the prediction error plots, 68% of all the times a 
certain prediction is made, the outcome will actually lie about a 
predicted mean, within a symmetric band having the half width indicated 
on the plot. 

Note that the regular rising and leveling trand of the prediction 
error plots is a sinusoidal form superimposed on a smoothly increasing 
curve. The sinusoid has a period of about 4 or 5 seconds. This is 
about one half the period of the ship motion. The same phenomenon 
also occurs for a single second order damped system and in the analysis 



141 


76265AW025 


2 


> 

t— 

Z 

< 


u 

Z 

D 


Z 

o 


O 


Q_ 



Prediction Errors for Vertical Ship Deck Position 
Figure 4. 3. 2. 2 
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Prediction Errors for Vertical Ship Deck Velocity 
Figure 4. 3. 2. 3 
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Predicted Correlation Coefficient for Ship Position and Velocity 
(Note: Irregularities are Due to Digital Calculation Errors) 

Figure 4. 3. 2. 4 



144 


of a sinusoidal wave with random amplitude A. The sinusoidal wave co- 
variance result is derived here. Let the sinusoidal wave be given as: 


y(t) = A sin (cut) (4.3.2.25) 

where the amplitude A is zero mean gaussian random variable, with 
variance a 2 . Then: 


E{ y(t) } = sin (cot) E {A} = 0 (4.3.2.26) 

Cov{ y(t) } = E{ [y(t) - E{ y(t)> ] 2 } = E{ y(t) 2 } 


= sin 2 (cot) E {A 2 } = ff 2 sin 2 (cot) 


s 


- a 2 (l - cos 2cot) 


(4.3.2.27) 


The squaring of the sinusoidal wave yields a frequency component in 
the covariance twice that of the original wave form. This corresponds 
to the observations made above for the ship motion. 

The ship vertical position and velocity correlation is shown in 
Figure (4. 3. 2. 4). The plot shows that position and velocity become 
less and less correlated in time. This means that knowing position 
or velocity will not generally help in determining the other. This is 
an important reason for rejecting landing control strategies that simply 
command hover at some fixed altitude until the ship rises high enough 
to contact the aircraft. Because ship position and velocity are un- 
correlated in steady state, the conditional distribution of ship 
velocities given ship height is the same at any height. Therefore, no 
altitude, even very high ones, will reduce the variation of ship 



145 


velocities at impact. Note that a simple second order damped oscillator 
also has this uncorrelated x, x property. For example, write the second 
order system as 

x = Ax + Bp (4.3.2.28) 


T 

where x » [x,x]' , and p is a zero mean gaussian white noise with inten- 
sity 1.0. Let the matrices A and B have the form: 


A = 


0 


- 0 ) 


2 


1 


-25u> 


(4.3.2.29) 


B = 


0 

< 1/2 


(4.3.2.30) 


Then the steady state solution of the state covariance equation 
(equation (4.3.2.22)) is given by: 


P 


q/4&> 3 0 

0 q/4£c») 


(4.3.2.31) 


where the two zero elements indicate that x and x are uncorrelated, 
as claimed. 


4.3.3 Nominal Trajectory, Design 

The purpose of the nominal trajectory is to guide the aircraft 
down to a predicted landing position using very little vertical control 
actuator effort. To accomplish this, a linear quadratic control prob- 
lem with a finite end time is set up and solved. The gains used to 
calculate a nominal trajectory during a landing are only dependent 
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on predicted time-to-go' and are stored as a function of that time. With 
the gains and nominal trajectory dynamics model available, only the ini- 
tial and final points and times are required to compute a nominal 
trajectory and its controls. 

The general form of the finite end time controller problem and 
solution are given as references (29) and (23). They are as follows: 
Find a feedback control law u(t) = ji(x(t)) to minimize 
the quadratic cost functional 


J(u) 




+ J F z T (t)Q(t)z(t) 


+ u (t)R(t)u(t)dt 


subject to 

the constraints of 

(4. 3.3.1) 

dynamics 

x(t) = Ax(t) + Bju(t) 

(4. 3.3.2) 

response 

Z(t) = Cx(t) 

(4. 3. 3. 3) 

initial conditions 2 £(tj) = x^ 

(4. 3. 3. 4) 

desired final conditions £(tp) =* 0 

(4. 3. 3.5) 

where (A,B) 

is stabilizable and (A,C) is detectable. 

F F = P Fi" 


is the quadratic cost associated with the final response y^t^), 

T 

Q(t) = Q (t) > 0 is the per unit time quadratic cost associated 

T 

with the response £(t), and R(t) = R (t) > 0 is the per unit time 
quadratic cost associated with the control .u(t) . 

The form of the feedback control _uQi(t)) is given as follows: 

_u(x(t)) = G(t)x(t) = -R ^(t)B T P(t)x(t) (4. 3.3.6) 

where P is the unique symmetric positive semidefinite solution of 
the matrix differential Riccati equation 



147 


P(t) = P(t)A + A T P(t) - P(t)BR“ 1 (t)B T P(t) +. Q(t) (4. 3. 3. 7) 

Note that the Rlccati equation is solved backwards in time from the 
end point, so the solution, and thus the gains can be stored as a 
function of " time-to-go" to the end. 

Using the assumption that aircraft vertical dynamics can be de- 
coupled from all other-degrees .of freedom, the -nominal, trajectory compu- 
tations are reduced to those for the third order vertical dynamics, 
which are given in Table (4. 3. 3.1). The states are deviation from the 
desired vertical end point x^, vertical velocity x^, and vertical 
acceleration x^. By the choice of the output matrix Cjj, the states 
are the outputs. The choice of control cost weighting R(t) is made 
unity as a reference level against which the output weighting Q(t) is 
compared in the cost functional. Q(t) has only one nonzero element 
weights the position. ■ has the following form: 

q N0M (t) = C0NST - ex P{ " - V 2/ a 2 } (4. 3. 3.8) 

The characteristic shape of ^Q^t) is shown in Figure (4. 3. 3.1). This 
form weights the position more heavily near the end of its trajectory 
than near the beginning, so that the feedback control solution will 
cause the deviation from the desired end position to approach zero 
smoothly. Note, can he written as a function of time~to^go, t^g, 

by setting t to t_ minus t__. q vn „(t) has three parameters. CONST is 
the scale factor indicating the importance of position error compared 
to the control effort, a is the parameter determining the width of the 


with terminal condition P(t^) 


NOMINAL TRAJECTORY DYNAMICS 


-N^ C) = + 

^(t) = c N % (t) = % (t) 


0 1.0 0 

0 0 32.2 

0 0 0 

0 

0 

1.0 

1.0 0 0 

0 1.0 0 

0 0 1.0 


Table 4. 3. 3.1 
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of the weighting in time and so indicates how soon the position error 
becomes important, and t^ is the time at which the error is weighted 
most heavily. The effect of making t^ < t^, is to cause' the velocity of 
the aircraft to be small as it reaches t^, for although position is 
weighted mostly at t , it is also weighted beyond indicating that 
the position error should remain near zero after t^. Since this weight- 
ing (t^ < tj.) causes the aircraft to reach the desired position and 
approximately zero velocity at t^, the gains obtained from solving this 
problem are shifted, so that those at time t^ are considered to be the 
final gains of a problem ending at t^. The gains from t^ to t^, are 
unused. The option of making t^ < t^ is referred to as "shifting the 
nominal trajectory gains'!. 

Two nominal trajectories are compared in Figures (4. 3. 3.2) and 
(4. 3. 3. 3). In these figures the nominal trajectory parameters are: 

CONST = 0.1 feet -2 

a = 1.0 second 

tp - t^ = 5.0 seconds 


Figure (4. 3.3.2) has a trajectory length of 6.0 seconds, while Figure 
(4. 3. 3. 3) has a trajectory length of 2.0 seconds. The short trajectory 
is shown to indicate how position errors 2.0 seconds before landing 
are reduced only partially, while in the longer one they were almost 
completely eliminated. 

The values of CONST, O and (t^, - t^) given above provided an 
initial set of nominal trajectory parameters for the design evaluation. 
They were selected because they produce smooth trajectories over time 
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Nominal trajectory for parameters t^, = 6.0, t = 11.0, 

a - 1.0 a * 


CONST = 0.1, 


Figure 4. 3. 3. 2 
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Nominal Trajectory for Parameters t^ = 2.0 t = 7.0 CONST = 0.1 a = 1.0 
Figure 4. 3. 3. 3 
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intervals on the order of 6.0 seconds. 

4.3.4 Trajectory Follower 

This section discusses the purpose of the trajectory follower, and 
the design used. 

The trajectory following controller has the objective of holding 
the aircraft to the commanded nominal trajectory in spite of any dis- 
turbances it may encounter. The corrections it makes to the aircraft 
trajectory are of higher frequency than those caused by the nominal 
commands. For this reason, the design of the follower is a steady 
state linear quadratic regulator with good tracking performance. Besides 
tracking the commanded vertical trajectory, the follower may also be 
used simultaneously to track ship orientation, lateral position, and 
fore-aft position, provided the necessary ship states are available. 

The control problem solved to get the feedback gains for the trajectory 
follower is identical to the one solved for the hover controller, with 
the exception that instead of weighting aircraft to ship vertical 

separation z , aircraft position alone was weighted by q z< The 

A 

weighting ratio q to q is the same as used previously for q Q 

9 A-S Z a 8 a-s 

to q . Based on the experience gained with the hover controller, 

Z A-S 

the quadratic weights were chosen as: 

q fl * io.o 

0 A-S 

q = 0.01 

Z A 

The gains obtained from these weightings are used to feed back vertical 
trajectory errors and aircraft to ship pitch orientation deviations to 
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the aircraft control inputs. 

4.3.5 End Point Attainability 

End point attainability refers to the ability of the controller to 
bring the aircraft to a given predicted end point using no more than a 
specified control authority level. Because of the random nature of the 
ship motions, not all deck motion crests are equally suitable for land- 
ing attempts. If a predicted crest is too low or occurs too soon, the 
nominal trajectory required to reach it may call for excessive control 
effort. Furthermore, as shown in Figure (4. 3. 3. 3), the nominal trajec- 
tory in some cases may not even reach the predicted end point. In 
other cases it may guide the aircraft to a high velocity impact with 
the ship. To prevent these situations, all nominal trajectories in- 
cluding ones computed as updates during a landing in progress should 
be checked for end point attainability. If control authority, miss 
distance, or impact velocity are unacceptable, a landing should not 
be initiated (a landing in progress should be aborted) . 

A crude example of such an attainability check is implemented in 
the present design. Each initial peak prediction is screened to deter- 
mine the control effort a bang-bang controller would use to move a 
point mass from the current aircraft altitude to the predicted height 
at the predicted time. Any end point requiring more than a specified 
control effort is skipped in favor of repeating the end point search 
at a later time. 

This simple end point attainability check is implemented with the 
window shown in Figure (4, 3. 5.1). Any point initially predicted in the 
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unshaded region is accepted and any outside the region is rejected. This 
screening is not repeated once a landing is underway. 

It should be noted that this screening does not mean the nominal 
trajectory computation will always use less than the specified actuator 
level, since the nominal trajectory minimizes a different cost functional • 
than the bang-bang controller trajectory (which minimizes time to the 
end given a maximum available control level) . However, the limit is a 
useful guide. 

This crude attainability check was found to be helpful during the 
design evaluation, but occasionally it proved inadequate. The algorithm 
checks only the initial prediction of the end point, and not the sub- 
sequent updates, which may be unattainable with the specified control 
actuator authority. A practical design should include a more complete 
checking scheme as originally described in this section. 

4.3.6 Update Determination 

If ship motion were deterministic, the state at the initial time 
would be sufficient to predict the end point perfectly, and thus the 
trajectory to reach it. However, this is not the case. Updates must 
be made to compensate for deviations from the predicted mean ship motion. 
The rate at which updates are made is a design decision within bounds 
imposed by two considerations. The upper rate is bounded by the 
frequency content of ship motion, while the lower rate is bounded by 
motion uncertainty. 
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First consider the upper bound on the update rate. This is deter- 
mined by the highest frequency components of ship motion. The Nyquist 
sampling frequency theorem, given in reference ( 3Q ) , states that the 
minimum sampling frequency to retain all the information contained in 
a band-limited signal is twice the maximum signal frequency. A sampling 
frequency higher than this is theoretically unnecessary. However, in 
practical implementations, where signals are not strictly band limited, 
a factor of ten increase beyond the minimum is often necessary. Since 
the dominant ship motion period is about 9.6 seconds, its frequency is 
0.65 radians /second. The corresponding Nyquist sampling frequency is 
1.3 radians/second, and with a factor of ten increase, a reasonable 
sampling frequency is 13 radians /second or once every 0.5 second. In 
the design evaluation both 0.5 and 0.2 second sampling periods are 
tested. 

The second consideration follows. The nominal landing trajectory 
is flat near the end point so any corrections made close to the terminal 
time can cause the aircraft position to track ship motion. If these 
corrections are made at a high rate, the landing controller mimics the 
performance of the hover controller. This performance has already been 
determined however, and found to be unsatisfactory since it uses high 
levels of control actuator effort. To avoid this tracking behavior, 
end point and nominal trajectory updates are suspended a certain time 
interval before the expected landing. This time is chosen as the 
largest possible, that does not incur a large ship motion prediction 
uncertainty. From the analysis of Section 4.3.2 on ship motion pre- 
diction, a value of 2.0 seconds was chosen. The value is somewhat 
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arbitrary, but it keeps the estimated mean ship velocity at touchdown 
near -0.8 feet/second as- listed in Table 4.3. 2.1. 

4.3.7 Overall Design 

The concept of the landing controller, with its main functional 
components has been described. In this section, the components are 
assembled in the overall desigi}. The type of evaluation chosen for 
the system performance is discussed also. 

A block diagram of the overall landing control loop is shown in 
Figure (4. 3. 7.1). It has three major blocks: one which generates the 

nominal trajectory x^ QM (t) and control one which computes 

the trajectory following feedback control 6u(t); and one which repre- 
sents the aircraft, ship and wave system. 

The nominal trajectory generator has four functions.. The first is 

t 

to determine when endpoint predictions and nominal trajectory updates 
are to be made. The second is to predict ship motion and identify its 
peaks. The next is to check end point attainability, and the last is 
to compute the nominal trajectory and control to the predicted end 
point. The nominal trajectory x^^Ct) and control u^^Ct) are the 
generator outputs, while the actual system state _x(t) at the update 
times is the generator input. 

The combined aircraft, ship, and wave system has two separate 
inputs. The wave receives disturbances p(t), which drive the wave 
dynamics, and hence produce stochastic ship motion. The aircraft 
receives commands _u(t), which are the sum of the nominal trajectory 
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Landing Controller Block Diagram 
Figure 4. 3. 7.1 
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and nominal trajectory tracking commands and <S_u(t) respectively. 

(Recall from Chapter Two on modelling, that wind disturbances acting on 
the aircraft have not been included to date.) The state vector jx(t) and 
the output vector _y(t) produced by the combined system contain aircraft, 
ship, and wave quantities. 

Finally, the trajectory follower computes the tracking command 
5_u(t), from the trajectory deviation 5x(t), which is the difference 
between the actual state :c(t) and the desired state x^Q^(t) . 

A sketch of the aircraft and ship tracks during a landing is drawn 
in Figure (4. 3. 7. 2). For clarity, the update rate in the sketch is 
only once every 2.0 seconds and the aircraft trajectory deviations are 
exaggerated. 

The evaluation of the overall controller design is discussed in the 
next section; however, before proceeding, an explanation of the form of 
the evaluation is appropriate. 

The system with the controller is nonlinear. In fact, both the 
dynamics and the output are nonlinear. First look at the nominal 
trajectory dynamics, 

= A N%OM ( - t) + (4. 3. 7.1) 

The nominal control u N0M (t) is 

“NOM^ = K tp^ %0M^ (4. 3. 7.2) 

so the closed loop dynamics are 


W C) = (A N + B H K tp (t)> W C) 


(4. 3. 7. 3) 
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where the feedback gain matrix K t p(t) is subscripted to denote that 

it is determined as a function of the time tp of the predicted ship 

motion crest, tp is in turn a nonlinear function of the ship and wave 

initial conditions. Thus, the true trajectory of the aircraft is a 

nonlinear fraction of the initial conditions. Secondly, look at the 

output. This is a vector composed of a linear combination of the. system 

state at the actual time of the touchdown t. : 

A 


H x(t A > 


(4. 3. 7. 4) 


t^ is a function of the ship-wave initial condition and its driving 
noise from the initial time to the landing time. Therefore, z. is also 
nonlinear in the initial condition since it depends on j£(t), which is 
nonlinear in x gilir ( £ j) (through K in the nominal trajectory) and on 
t^, which is also dependent on x ^ T p (t T ) . Because of these nonlineari- 
ties, Monte Carlo simulations are chosen for the evaluation. 


4.4 Design Evaluation 


As explained above, the landing controller was evaluated by 
direct Monte Carlo simulations. The following sections discuss its 
performance under various parameter modifications , and compare its 
performance with the hover controller of Chapter Three. The- last 
section summarizes the evaluation and comparison results. 
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4.4.1 Evaluation Parameters 

The effects of five parameters on the landing performance were 
investigated. These were aircraft initial hover height, endpoint and 
trajectory update rate, nominal trajectory gain shift, update suspension 
time, and control actuator nominal authority. Table (4.4.1!.l) contains 
the numerical results of the Monte Carlo simulations. 

The results for aircraft initial hover heights of 13 and 8 feet 
above the mean deck height are shown at the top of the table. In 
this comparison, the update suspension is not implemented. For these 
heights, the average impact velocity decreases with decreasing height. 

In addition, the actuator level required decreases by a factor of two. 

The decrease in actuator effort is a significant amount and is 
consistent with the physics of the problem: to displace a mass a long 
distance in a fixed time, with zero initial and final velocities, re- 
quires more effort than to displace it a shorter distance. 

The second comparison in Table (4. 4. 1.1) shows that an increased • 
rate of prediction and trajectory updates has the effect of decreasing 
impact velocity. This is because updates occur nearer the end time 
and thus less prediction errors buildup before touchdown. With only 
an initial prediction and no updates, the impacts have average velocity 
near 3 feet/second. When updates every 2 seconds are implemented, 
this velocity is reduced to under 2 feet/second. More frequent updates 
continue to reduce the impact velocity. At a rate of one update every 
0.2 second, the velocity is 1.6 feet/second. Again, this evaluation 
does not include the update suspension feature. 


Evaluation: The Effects of Five Parameters on Landing Performance 
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Nominal trajectory gain shifts of 5 and zero seconds are compared 
next, with and without update suspension. The results show that without 
a gain shift, the average impact velocity is lowest but the variation of 
impact velocity about the mean is large. This tradeoff between average 
performance and variation about the average is independent of update 
suspension. 

Update suspension before the predicted end time yields higher impact 
velocities than no update suspension does, whether or not nominal tra- 
jectory gains are shifted. However, in this comparison shifted gains 
do reduce the required maximum actuator effort. 

The final comparison in Table (4. 4. 1.1) shows that control actuator 
nominal authority limitation reduces actuator excursions and the average 
aircraft-ship impact velocity. The actuator excursion reduction is 
especially significant. Further reductions appear possible if the 
nominal trajectory control sequence is assessed at each update rather 
than at the initial time only. 

The set of landing controller parameters, which yielded the lowest 
impact velocities during the evaluation simulations, are as follows: 
initial height -13.0 feet 

update period 0.2 seconds 

nominal trajectory gain shift 0.0 seconds 

update suspension interval before the predicted time 0.0 seconds 
control actuator nominal authority 3.0% of nominal 

For this set of parameters, the average impact velocity was about 1.0 
feet/second with a standard deviation of about 0.7 feet/second. These 
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statistics were computed from eleven Monte Carlo simulations. During 
these eleven simulations, the maximum negative actuator effort was 
-6.2 percent of nominal. (Note, that the maximum negative value is 
of interest since this represents a total thrust exceeding the hover 
level. Positive actuator effort represents a reduction in overall 
thrust.) 

Figure (4.4. 1.1) shows time histories of a typical landing using 
the above controller parameter values. The results obtained are the 
following: 

relative aircraft-ship impact velocity 1.3 feet/second 

vertical aircraft velocity -0.3 feet/second 

vertical ship velocity -1.6 feet/second 

vertical ship position -5.9 feet 

maximum actuator effort employed 2.9% of nominal 

An examination of the timing and height prediction errors made 
on the final update before touchdown reveals that landings occur earlier 
and at more elevated positions than predicted. Figure (4. 4. 1.2) 
displays the experimental distribution of errors obtained from eleven 
simulations using the same controller as in the sample landing above. 

The figure shows that landings occurred between 0.2 and 0.7 seconds 
earlier and between 0.0 and 0.8 feet higher than the final predictions 
made before the touchdowns . These timing errors agree qualitatively 
with those calculated in the analysis of Section (4.3.2), although they 
are significantly larger. It is likely that this discrepancy is attri- 
butable to the simplifying assumptions and approximations made in 
the analysis . 
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4.4.2 Comparison with the Hover Controller 

The results for the landing controller indicate better performance 
than the hover controller. Using the best set of parameters as indicated 
in the previous section, the landing algorithm yields an average touch- 
down velocity of about 1.0 feet/second, whereas the hover controller, 
which uses the same control actuator effort, yields an average touch- 
down velocity of about 1.8 feet/second. Alternatively, to obtain the 

\ 

same average impact velocity as the landing controller, the hover con- 
troller must use about twice as much control authority as the landing 
controller. 

This comparison does not take into account the practical consi- 
derations required for implementation. These additional considerations 
also favor the landing controller since it requires less ship motion 
data and actually computes the controls necessary to bring the air- 
craft from stationary hover, along a nominal trajectory, to a shipboard 
landing. 

The comparison is presented graphically in Figure (4. 4. 2.1). The 
figure shows the distribution of impact velocity and extreme negative 
excursions of the vertical control actuator (in terms of thrust to 
weight ratio) obtained for the simulations. For comparison with the 
hover controller, the figure includes the mean and one rms bands for 
the hover controller impact velocity. The computation of the mean and 
standard deviation for the hover controller is described below. 

Aircraft-ship touchdown only occurs when the relative separation 
z^_g(t) is zero and the relative velocity z^_^ (t) is positive. For 
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the hover controller these quantities are stationary, independent 

gaussian random processes with zero means, and standard deviations 

a and a. , respectively. (In figure (3.4. 3.2), the three-sigma 
Z A-S Z A-S 

values are plotted as functions of thrust to weight.) The expected 
touchdown velocity is given by: 

E {z A-S (t) I Z A-S (t) = °» Vs (t) - 0} (4. 4. 2.1) 

but the distributions of z^_^ (t) and z A _<,(t) are independent of time, 
and also statistically independent of each other, so the expectation 
reduces to : 

E{ z A _ s | z A _ s >_ 0 } (4. 4. 2. 2) 

The probability density function for the velocity v =* z A _g is: 

P (V) =- (l//2ir) (1/a ) exp{ --j(V/a v ) 2 } (4. 4. 2. 3) 

for — 00 < V <°°. Hence, the conditional density is: 

P v | v > 0 (v|v > 0) - (/2Af)(l/a v ) exp{- ~r (V/a v ) 2 } (4. 4. 2. 4) 

for 0 _< V < 00 . The resulting mean and standard deviation of the impact 
velocity are: 

E{ Z A-S I Z A-S — 0 1 = (A/tt ) a* (4. 4. 2. 5) 

A-s 

( Var{ z A _ s | z A _ s >. 0 } ) ly ^ 2 » (A -2/ir)cr» (4. 4. 2. 6) 

A-S 

Using the conversion between a* and thrust to weight available from 

2 A-S 

Figure (3. 4. 3. 2), this mean and standard deviation are plotted in Figure 
(4. 4. 2.1) for comparison with the landing controller simulation results. 
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4.4.3 Conclusions 

The goal of the landing controller was to decrease the thrust to 
weight ratio required to attain low aircraft-ship relative touchdown 
velocities compared to the value required by the tracking hover con- 
troller. This goal was achieved in Monte Carlo simulations of the 
controller. 

The landing controller achieved touchdown velocities near one 
foot/second. These were obtained with a 0.2 second sampling rate and 
a thrust to weight control authority approximately half the magnitude 
required by the hover controller yielding the same average impact 
velocity. The extreme values of thrust to weight were 5 or 6 percent 
above the stationary hover value of 1.0. High initial hover positions 
required higher thrust to weight ratios. The nominal trajectory gain 
shift did not yield the expected softer landings that the earlier 
nominal trajectory analysis predict. Similarly, update suspension 
before touchdown did not reduce impact velocity or thrust to weight 
requirements. Finally, control acutator nominal authority limitation 
proved to be beneficial by significantly reducing the control levels 
used. 

These observations indicate that the best performing landing con- 
troller is one using, unshifted nominal trajectory gains, no update 
suspension, minimum safe initial hover height and update rates of 5 
per second or higher. 
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The biases between the actual and predicted touchdown itmes and 
heights that were found analytically and also observed in the simula- 
tions, result in significant touchdown velocities . 

Unlike the hover controller, the landing controller performance 
is not necessarily the "best" possible. Other landing algorithms may 
yield better performance, either by using lower control power or ob- 
taining lower touchdown velocities. This is because the landing con- 
troller algorithm is only an approximate solution of its defining 
optimization problem in Section (4.2.1). On the other hand, physical 
implementation of the algorithm probably would not perform as well as 
the simulation results indicate, since the simulations were based on 
idealistic modelling assumptions. In this sense, the landing control- 
ler results presented in this chapter can still be viewed as lower 
performance limits fo'f VTOL systems. 
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CHAPTER 5 
CONCLUSIONS 


5.1 Introduction 

The pr ima ry objective of this research has been to design a ship- 
board landing controller for a hovering VTOL aircraft. The major prob- 
lems faced in such a landing, especially aboard a small ship in heavy 
seas, are the limited control authority and the severe ship motions. 

When these difficulties are overcome, VTOL operations from small ships 
will be considerably closer to realization than they are currently. The 
results of this research are basic performance characteristics and 
requirements for VTOL controllers. They are not intended as detailed 
engineering designs . 

5. 2 Achievements 

Aircraft dynamics, wave motions, and ship motion were modelled. 

The resulting overall model was then used in the design and analysis of 
the hover and landing controllers . The hover controller demonstrates 
that good shipboard landing performance is possible. The landing con- 
troller exhibits similar performance levels, but with reduced control 
power requirements . 

The aircraft considered is a representative lift/cruise fan VTOL 
modelled near stationary hover. Each of the six degrees of freedom 
associated with the rigid vehicle motion are controlled by a separate 
actuator with assumed first order dynamics. Aerodynamics were not 
modelled. However, the model is structured to accommodate changes and 
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additions such as aerodynamics for future studies. Sea state five ocean 
wave conditions were modelled as shaped white noise with a theoretical 
Neumann spectrum, modified to account for ship speed and heading. The 
ship itself is a small destroyer of 400 foot length. Its pitch and heave 
and roll motions are driven by two identical, but independent wave models 
Each of the three motion responses was modelled as a second order system, 
with parameters estimated from available frequency response data of 
similar ships. 

Two VTOL control concepts were designed and analyzed. The first is 
a "hover controller", which commands the vehicle to track ship landing 
pad motions with small control actuator thrust levels. A steady state 
linear quadratic regulator design was used for this controller. The 
tracking errors, control effort, and required sensor bandwidth were 
determined analytically. For 3 rms tracking errors on the order of one 
foot, thrust to weight ratios near 1.2 are required by this controller. 
These results represent ultimate (lower) performance bounds, since not 
all the dynamics, disturbances, and measurement constraints are accounted 
for. The second design is a "landing controller". Its goals are to 
land the vehicle onboard the ship with a low touchdown velocity and small 
control thrust levels. The "landing controller" commands the aircraft 
from stationary hover along a smooth, low control effort trajectory, to 
a touchdown on a predicted crest of ship motion. Performance analysis 
was by Monte Carlo simulations since the overall system is nonlinear; 
and hence, incompatible with linear covariance analyses. The effect 
of five algorithm parameters were studied. With "optimized" parameters. 
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the landing controller uses about a half of the control authority 
required by the hover controller, which has comparable touchdown velo- 
cities. For average touchdown velocities of one foot per second, the 
thrust level requirements are about six percent of the nominal hover 
amount. These landing controller results are subject to the same per- 
formance bound interpretation as the hover controller results since the 
same aircraft-wave-ship system model assumptions were used. 

i 

5.3 Further Research 

As this work stands, several open questions remain. These concern 
the influence of the following: 

(1) the inclusion of full longitudinal and lateral dynamics, 

(2) noisy and incomplete sensor measurements, 

(3) human piloting considerations, 

(4) ship and ocean conditions, 

(5) aerodynamic disturbances, and 

(6) ground effects. 

The first of these is of considerable interest in the light of 
the low levels of lateral control that have been displayed by VTOL 
aircraft, which limits their maneuverability ( 31 ) . The influence of 

sensor measurements should be examined, since it is improbable that all 
system states will be available for feedback. Particularly, it is im- 
probable that wave states can be measured. Therefore, the problem of 
predicting ship motion based only on practical measurements is an im- 
portant research question. Pilot dynamics representing human tracking 
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of ship motion or of a commanded trajectory can be incorporated into 
the overall model. A comparison of the subsequent results with those 
obtained here would indicate the performance degradation associated with 
inserting, a pilot into the loops of the automatic hover and landing con- 
trollers. If. was found in the hover controller performance analysis, 
that ship motion has a strong effect on performance. Hence a parametric 
study of controller performance versus different ship types, speeds, 
headings, and sea states could yield useful guidelines for controller 
design and general operating procedures for VTOLs at sea. Aerodynamic 
forces and moments on the aircraft are unlikely to be negligible at the 
high wind speeds crossing the ship deck. As a result, their effects 
should be evaluated. VTOL vehicles are also susceptible to ground 
effects such as fountain and suck-down phenomena. In the past, these 
have not been reliably predicted (4 ) . However, by making a simple 

model of the force and moment magnitudes they exert on the aircraft as 
a function of, say, aircraft to landing deck separation and orientation, 
helpful parametric information relating controller' requirements to 
ground effect magnitude may be obtained. A study of this effect is 
strongly recommended, since it is critical to VTOL aircraft operating 
in ground proximity. 

This work has laid a foundation for these studies. A versatile 
model and promising modern control design and analysis approach have 
been developed for shipboard landing of VTOL aircraft. From this point, 
more accurate and more complete models may be used with the same design 
and analysis techniques, or better performing control algorithms may be 
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developed, analyzed and compared with the current model and results 
established in this report. 
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APPENDIX 2. 4. 3.1 


MATRICES USED IN THE AIRCRAFT MODELLING 
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APPENDIX 2. 4. 3. 2 


DERIVATION OF RAM DRAG FORCE AND MOMENT EXPRESSIONS 


Ram drag forces and moments are those caused by the change in linear 

♦ 

momentum given to the air flow at fan and engine inlets because the air- 
craft is translating and rotating. (The momentum change due to acceleration 
of air in the fans that generate nominal control thrusts is. accounted for 


in the rigid body dynamical equations.) Below, a short derivation of 
the ram drag force and moment equation is given. Note that these equa- 
tions are for the forces and moments f ^ ^ and acting on the air. 

The complementary ones f^ and m^ acting on the aircraft are in the 
opposite directions. The final equations are summarized first. 
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nu is the airmass flow through fan or engine i and 


In these equations 
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and is the position of that fan or engine inlet with respect to the 
aircraft center of mass. 

To begin the ram drag force derivation, Newton's Law (or the prin- 
ciple of linear momentum conservation) is written in the earth frame 
for a mass m. : 

i 

hi * 4- hi. “ 4- (, lV ‘ + "Ai (A2.4.3.2.3) 

dt dt 


For a rigid body, m^ = 0; however, the engines and fans are not rigid 
bodies since there is mass flow m through them. (The rotating parts are 
considered as contributing internal angular momentum in the aircraft 
moment equation.) The m^v^ term is nonzero and is refered to as the 
ram drag force. 


■^REFi 


m, v_. 
i— El. 


(A2.4.3.2.4) 


This is the expression used for the ith fan or engine inlet. In body 
coordinates, were z_^ is the location of the ith inlet, the force is: 


^SBFi 




+ ^ x r.) 


(A2.4.3.2.5) 


Summing over all the inlets gives equation (A2.4.3.2.1) again: 


■^RBF 




(A2.4.3.2.6) 


The conservation of angular momentum can also take into account the air 
mass flows. A torque or moment m. on a small mass m. at distance z, 
from the center of rotation (assumed to be the center of mass here) is 
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SaBF ■ “r^*2b + (A2.4.3.2.10) 

which is the last of the two equations describing ram drag effects, 
equation (A2.4.3.2.2) . 

Although, in general, the nu are time varying, they are considered 
constant in the aircraft model. 
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